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ABSTRACT 


The  ability  to  make  computer  images  more  realistic  is  becoming  more 
important  as  the  hardware  for  producing  such  images  is  becoming  less  expensive 
and  hence  more  available.  The  key  to  producing  realistic  images  lies  in  the 
algorithms  that  can  take  full  advantage  pf  the  hardware  to  produce  them.  In  this 
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study,  we  look  at  a  prototype  of  i*ray  tracer^  as  presented  in  [Ref.  l).  Ray  tracing, 
in  combination  with  a  global  illumination  model,  currently  provides  the  most 
realistic  images  that  can  be  generated  on  general  purpose  computing  hardware. 
The  prototype  was  succesfully  implemented  on  an  IBM  AT  clone. 
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From  the  beginning  of  recorded  history  mankind  has  always  had  the  need  to 
create  pictures.  The  reasons  for  creating  these  pictures  ranges  from  the  aesthetic, 
pretty  pictures  are  nice  to  look  at,  to  the  functional,  pictures  can  be  an  excellent 
way  to  communicate  information.  As  mankind  progressed,  so  did  his  ability  to 
create  pictures,  although  the  techniques  used  to  create  pictures  basically  stayed 
the  same.  The  advent  of  computers  gave  man  yet  another  tool  with  which  to 
create  pictures. 

The  rapid  increase  in  technology  has  made  computer  graphics  a  rapidly 
growing  field.  For  the  first  time  since  man  started  drawing  pictures,  completely 
new  techniques  needed  to  be  developed.  In  computer  graphics  the  brush,  paint 
and  canvas  are  replaced  by  the  mouse,  algorithm  and  display.  However  even 
though  the  tools  have  changed,  the  same  problems  remain:  how  to  make  the 
picture  look  better,  be  it  either  more  pleasing  to  the  eye  or  to  get  the  information 
across  more  clearly. 

Two  of  rhe  most  common  and  difficult  problems  in  computer  graphics  are  'he 

hidden  surface  and  lighting  and  shading  problems.  A  large  number  of  solutions 
exist  to  both  of  these  problems.  Very  few  solutions  can  be  applied  to  both.  One 
such  solution  is  ray  tracing.  Ray  tracing  is  the  process  of  following  an  imaginary 


ray  from  a  viewpoint  through  a  pixel  on  a  screen  and  into  a  scene  to  determine  if 
it  intersects  any  objects  in  the  scene  and  then  calculating  the  intensity  of  the  pixel 
it  went  through  based  on  the  the  final  destination  of  the  ray.  As  in  most  cases 
where  one  solution  is  found  to  several  problems,  that  solution  is  seldom  the  best 
for  all  the  problems  it  is  applied  to  and  so  it  becomes  a  matter  of  trade-offs.  Such 
is  the  case  with  ray  tracing.  Amorg  the  hidden  surface  removal  techniques,  ray 
tracing  is  the  least  efficient  being  referred  to  as  a  brute  force  technique.  In 
contrast,  it  has  been  labeled  as  one  of  the  most  elegant  techniques  in  regards  to 
lighting  and  shading  [Ref.  2:  p.  137],  Because  of  this  latter  fact,  ray  tracing  has 
become  an  important  technique  in  computer  graphics.  Ever  since  the  idea  behind 
ray  tracing  was  suggested  by  Appel,  numerous  articles,  studies,  and 
implementations  have  been  done  on  it.  These  in  turn  have  spawned  fruther 
extensions  and  modifications.  [Ref.  1:  p.  296] 

A.  DEFINITION  AND  OVERVIEW 

The  idea  behind  ray  tracing  lies  in  the  theory  that  the  light  in  our 
environment  can  be  modeled  as  rays.  After  being  emitted  from  a  source,  the  rays 
are  then  reflected  and  refracted  through  a  scene.  Some  of  the  rays  eventually  find 
•  heir  vay  *o  *ne  eve  where  *he  scene  is  recreated  Figure  1.1a:.  These  iignt  rays 
are  emitted  from  light  sources,  such  as  the  sun.  An  infinite  number  of  light  rays 
exist,  but  only  a  small  percentage  of  them  are  received  by  us.  To  try  and  trace 
these  rays  from  the  source  is  computationally  expensive.  Appel  suggested  that 
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Tracing  Rays  Backwards  from  View  Position 


instead  of  tracing  the  rays  from  the  source  that  they  should  be  traced  backwards 
from  the  viewer,  thus  dealing  with  only  those  rays  that  actually  contribute  to  the 
scene  (Figure  1.1b).  [Ref.  1:  p.  296] 

The  basic  ray  tracing  algorithm  is  a  very  simple  one  and  not  difficult  to 
implement.  The  basic  algorithm  is  a  hidden  surface  algorithm.  All  hidden  surface 
algorithms  can  be  classified  based  on  the  coordinate  system  or  space  in  which  they 
are  implemented.  These  are  either  in  object  space  or  in  image  space.  The  ray 
tracing  algorithm  falls  under  the  category  of  image  space.  This  category  of 
algorithm  is  implemented  in  the  screen  coordinate  system  in  which  the  objects  are 
viewed.  Unfortunately,  the  calculations  are  performed  only  to  the  precision  of  the 
scene  representation,  which  generally  provides  poor  resolution.  The  image  space 
algorithms  work  by  comparing  every  object  in  the  scene  with  every  pixel.  Such  an 
algorithm  is  computationally  expensive.  Ray  tracing  algorithms  have  three  parts: 
a  viewpoint,  a  raster  screen,  and  a  set  of  objects  (Figure  1.2).  In  the  algorithm, 
the  viewpoint  is  along  the  positive  z  axis.  From  this  point,  a  ray  is  shot  into  the 
scene  through  the  center  of  every  pixel  on  the  raster.  Each  of  these  rays  is  then 
traced  and  compared  against  every  object  in  the  scene  to  determine  if  there  is  an 
intersection  with  any  of  them.  It  is  in  the  determination  of  a  possible  intersection 
poinr  "hat  i  *ay  ‘racer  -spends  anywhere  from  75  -o  95  percent  of  its  rime,  if 
there  is  an  intersection,  then  the  intensity  at  the  pixel  is  determined  using  the 
intersected  object’s  attributes  and  an  appropriate  illumination  model.  If  there  is 
no  intersection,  then  the  pixel  intensity  is  determined  by  the  background 


F igure  1 . 2 

Layout  of  Ray  Tracing  Scene  [Ref.  1:  p.297] 
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intensity.  This  procedure  is  then  repeated  for  every  pixel  on  the  raster.  When 
the  ray  tracer  is  used  as  a  hidden  surface  algorithm,  intersection  testing  stops 
after  the  first  intersection.  Extensions  to  the  basic  ray  tracing  algorithm,  that 
showed  its  usefulness  in  implementing  a  global  illumination  model,  were  originally 
implemented  by  Whitted  [Ref.  4]  and  Kay  [Ref.  5  and  6].  In  these  extensions  of 
ray  tracing,  additional  rays  are  calculated,  specifically  the  reflected  and  refracted 
rays,  and  then  tested  to  see  if  they  intersect  with  any  objects  in  the  scene.  This 
process  of  generating  new  rays  and  tracing  them  to  check  for  possible  intersections 
is  continued  until  the  rays  either  leave  the  scene  or  stack  space  is  exceeded.  In 
such  a  case,  the  remaining  rays  are  treated  as  if  they  had  left  the  scene.  This 
process,  illustrated  in  Figure  1.3a,  for  a  single  ray  with  intersections  is  easily 
represented  using  the  tree  structure  shown  in  Figure  1.3b.  Here  each  node  of  the 
tree  represents  a  ray  surface  intersection.  At  each  node,  at  least  one  and 
sometimes  two  subbranches  are  generated.  One  branch  of  each  of  the  reflected 
and  refracted  rays  is  generated  from  the  point.  [Ref.  1:  pp.  190-296] 


B.  ORGANIZATION 

This  study  is  broken  into  three  areas:  data  requirements,  ray  tracing 
me'houoiogy.  anti  ‘he  intensity  problem.  The  first  section  reviews  the  data  needed 
for  a  lighting  and  shading  modeler,  hereafter  referred  to  as  a  Tenderer,  of  which  a 
ray  tracer  is  an  integral  part.  The  second  section  reviews  the  actual  process  of 
tracing  a  ray  through  a  scene  to  be  rendered.  The  third  section  looks  briefly  at  the 
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illumination  problem  and  how  it  relates  to  the  ray  tracer.  The  concluding 
chapters  present  the  implementation,  and  known  limitations  of  the  model  along 
with  areas  of  future  research. 


II.  DATA  REQUIREMENTS 


A.  OVERVIEW  OF  THE  DATA  REQUIREMENTS 

The  importance  of  the  intersection  routines  in  the  ray  tracer  is  apparent  in 
the  fact  that  a  ray  tracer  spends  75  to  95  percent  of  its  time  determining 
intersections  [Ref.  1:  p.  297].  The  key  to  determining  intersections,  however,  lies 
in  large  part  on  the  data  used  to  describe  the  scene  that  is  being  rendered. 
Information  is  needed  not  only  to  describe  the  entire  scene  that  is  being  rendered 
but  more  importantly  to  describe  each  object  in  the  scene.  Scene  data  is  that 
information  needed  to  completely  describe  a  picture,  i.e.,  the  number,  kind,  shape, 
and  color  of  any  objects  in  the  picture  along  with  the  background  intensity  and 
light  source  information.  This  information  must  be  properly  ordered  and  broken 
down.  Falby  [Ref.  3]  suggested  that  a  scene  be  broken  into  three  categories: 
object,  view,  and  light.  Each  of  these  areas  is  examined  below  in  the  context  of  a 
ray  tracing  algorithm. 

1.  Object  Data 

The  data  pertaining  +o  each  object  in  ‘he  scene  can  be  grouped  into  +wo 
categories;  polygon  and  oounding  volume.  The  reason  lor  this  oreakdown  ;s 
twofold.  First,  each  object  in  the  scene  is  composed  of  polygons.  They  are  the 
basic  building  blocks  of  the  scene.  Second,  in  order  to  reduce  the  number  of 


intersection  checks,  it  is  necessary  to  set  up  a  boundary  around  each  object  so 
that  the  ray  tracer  only  performs  the  intersection  checks  in  the  regions  that 
actually  contain  an  object.  Such  a  boundary  is  called  a  bounding  volume.  We 
examine  the  polygon  data  first, 
a.  Polygon  Data 

Since  the  main  focus  of  the  ray  tracing  algorithm  lies  in  determining 
the  intersections  between  the  rays  shot  into  a  scene  and  the  objects  that  make  up 
a  scene,  and  since  each  object  is  comprised  of  polygons,  the  problem  is  really  one 
of  determining  the  intersection  points  between  the  rays  and  the  polygons.  From 
the  fundamentals  of  vector  calculus,  it  is  known  that  in  order  to  determine  the 
intersection  between  a  ray  and  a  polygon  only  the  vertices  of  the  polygon  are 
needed  as  well  as  the  direction  of  the  ray  and  a  point  on  the  ray.  Since  the  object 
is  to  be  constructed  of  polygons,  its  vertices  are  known.  Therefore,  it  is  only 
necessary  to  ensure  that  these  points  are  stored  in  some  manner,  such  as  a  record, 
so  as  to  be  accessible  to  the  ray  tracer.  In  order  to  determine  what  the  intensity 
of  the  pixel  is  through  which  the  ray  passes,  it  is  essential  that  the  characteristics 
of  the  object  whose  polygon  was  intersected  be  available.  Since  an  object  is  made 
of  polygons,  they  inherit  the  characteristics  of  the  object.  These  characteristics 
•tiso  need  ro  ’>e  readily  accession*  uid.  herefore,  iced  *o  be  stored  in  ->oine 
manner.  The  following  is  a  list  of  the  basic  object  characteristics  that  need  to  be 
available:  (1)  the  specular,  diffuse,  and  transmission  coefficients;  (2)  the  Phong 
specular  exponent;  and  (3)  the  index  of  refraction,  see  Table  2.1.  [Ref.  3:  p.  68] 


TABLE  2.1:  OBJECT  DATA 


FIELD  NAME 

VARIABLE  NAME 

VALUE 

Polygon  Vertices 

x,  y,  z 

real 

Diffuse  Coefficient 

Kd ,  . 

real 

:r0-l) 

Specular  Coefficient 

real 

Transmission  Coefficient 

KtrlS,6 

real 

(o-i) 

Unit  Surface  Normal 

z,  y,  z 

real 

ks» 

Phong  Specular  Exponent 

n 

Index  of  Refraction 

* 2 

real 

b.  Bounding  Volume  Data 

The  major  disadvantage  of  ray  tracing  is  that  it  takes  so  much  time. 
This  is  hard  to  avoid  since  it  is  so  computationally  expensive.  It  is  essential, 
therefore,  that  more  efficient  techniques  be  developed  to  assist  in  reducing  the 
number  of  calculations.  Several  techniques  already  exist  with  the  bounding 
volume  being  the  most  effective  [Ref.  1:  p.  298].  In  the  description  of  ray  tracing 
given  so  far,  it  has  been  stated  that  a  ray  is  checked  to  see  if  it  intersects  with  any 
object.  Upon  dissecting  this  statement  further,  a  better  understanding  of  the 
intersection  problem  can  be  realized.  Unless  some  optimization  is  done,  the  ray 
tracing  algorithm  is  forced  to  do  the  following.  Each  ray  must  be  checked  for  a 
possible  intersection  with  each  object.  Since  each  object  is  made  up  of  polygons, 
then  the  ray  must  be  checked  for  a  possible  intersection  with  each  polygon.  For  a 
■ornpiicared  object,  such  as  a  "eapor.  this  requires  a  large  number  of  checks  and 
must  be  done  for  each  object.  The  purpose  for  establishing  the  bounding  volume 
lies  in  two  facts.  The  first  is  that  generally  scenes  are  mostly  background  with 


just  a  few  objects,  hence  very  few  of  the  rays  actually  hit  anything.  Therefore, 


most  of  the  intersection  tests  done  are  a  waste  of  time.  Second,  a  ray  can  only  hit 
one  object  at  a  time.  To  have  it  process  through  the  entire  list  of  objects,  when 
intersections  with  most  of  them  can  be  eliminated,  is  needless.  A  bounding 
volume  is,  therefore,  a  method  of  enclosing  each  object  in  the  scene  in  a  simple 
containment  vessel,  which  in  effect  creates  a  boundary  around  the  object.  Once 
this  boundary  is  established,  the  number  of  overall  intersection  tests  can  be 
greatly  reduced,  as  in  the  example  of  a  teapot,  which  might  easily  have  over  a 
hundred  polygons.  If  it  is  surrounded  by  a  bounding  box  consisting  of  just  six 
polygons,  the  number  of  intersection  tests  can  be  significantly  reduced.  In  this 
situation,  instead  of  having  to  test  each  ray  against  each  polygon  of  the  object, 
only  those  rays  that  penetrate  the  bounding  volume  need  to  be  checked.  Thus 
the  bounding  volume  is  a  way  to  filter  out  unnecessary  intersection  tests  by 
limiting  the  tests  to  those  rays  that  are  most  likely  to  intersect  an  object. 

Just  as  the  use  of  a  bounding  volume  greatly  increases  the  efficiency 
of  the  ray  tracing  algorithm,  the  use  of  the  right  kind  of  a  bounding  volume  can 
improve  upon  that  even  more.  In  Rogers  [Ref.  l],  the  bounding  volumes 
suggested  are  a  bounding  box  and  a  bounding  sphere  (Figure  2.1),  each  of  which 
has  advantages  and  disadvantages.  The  bounding  sphere  is  much  easier  to 
implement  although  it  is  less  efficient  in  reducing  the  *arget  area  'hun  is  the 
bounding  box,  see  Figure  2.1.  The  bounding  box,  on  the  other  hand,  is 
computationally  expensive  to  implement.  The  data  needed  to  establish  a 
bounding  sphere  is  minimal.  It  only  requires  a  center  point  for  the  object  and  a 


radius  that  encompasses  every  point  of  the  object.  The  bounding  box,  on  the 
other  hand,  requires  far  more  data  in  that  the  polygons  that  make  up  the  box 
must  be  described. 


2.  View  Data 

In  rendering  any  scene,  certain  information  can  be  applied  to  the  scene 
as  a  whole.  This  information  is  grouped  together  to  form  the  view  data.  This 
data  consists  of  the  viewpoint  position,  a  constant  to  prevent  division  by  zero,  a 
refraction  index  for  the  global  medium,  the  ambient  light  intensity,  the 
background  color,  and  the  scene  dimensions,  see  Table  2.2.  [Ref.  3:  pp.  74-75] 

3.  Light  Data 

To  support  a  lighting  and  shading  model,  it  is  necessary  to  include 
certain  information  on  the  light  source  for  the  scene.  That  information  must 
include  the  position  of  the  light  source,  its  intensity,  its  type,  geometry,  and 
dimension,  see  Table  2.3. 


TABLE  2.2:  VIEW  DATA 


TABLE  2.3:  LIGHT  DATA 


FIELD  NAME 

VARIABLE  NAME 

VALUE 

Light  Position 

x,  y,  2 

real 

Intensity 

r 

real  (0-1) 

Type 

type 

enumerated  (point, distributed) 

Shape 

shape 

enumerated  (circular,  rectangular) 

Dimensions 

x,  y 

real 

B.  DATA  STRUCTURE  FOR  A  RAY  TRACER 

Falby  [Ref.  3]  suggested  a  data  structure  for  a  multi-illumination  model 
Tenderer.  That  data  structure,  with  minor  variations,  has  been  used  in  this  study. 
In  [Ref.  3]  a  complete  derivation  of  the  data  is  presented,  so  for  the  purposes  of 
this  study  only  a  brief  description  is  given  here.  Figure  2.2  illustrates  the  layout 
of  the  data  structure  as  used  in  this  study.  This  data  structure  essentially  consists 
of  arrays  of  records  layed  out  in  a  hierarchical  organization.  Starting  from  the 
highest  level  it  consists  of  the  following:  a  picture  record,  an  array  of  light  records, 
an  array  of  objects,  an  array  of  subobjects,  an  array  of  common  part  records,  an 
array  of  polygons,  and  three  arrays  for  the  vertices.  Each  of  these  is  now 
examined. 

1.  Picture 

Picture  is  a  single  record  which  contains  the  view  data  mentioned  earlier. 

2.  Lights 

The  lights  array  is  an  array  of  records,  with  one  record  for  each  light 


source  in  the  scene.  Each  record  contains  the  light  data  mentioned  above. 
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3.  Objects 


The  objects  array  is  an  array  of  records,  with  one  record  for  each  object 
in  the  scene.  In  this  study,  an  object  is  the  highest  order  item  in  a  scene.  Just  as 
the  scene  is  divided  up  into  objects,  each  with  its  own  bounding  volume,  so  is 
each  object  broken  down  into  subobjects,  each  with  its  own  bounding  volume. 

4.  Subobjects 

The  subobjects  array  is  an  array  of  records,  with  each  array  belonging  to 
one  object  record.  For  example  Figure  2.3a  shows  one  object,  a  barbell,  that  is 
divided  into  three  subobjects  which  are:  the  left  weight,  the  right  weight,  and  the 
bar.  The  record  layout  for  this  is  as  illustrated  in  Figure  2.3c.  A  subobject  is  the 
smallest  item  in  the  scene.  Each  object  has  at  least  one  subobject.  A  subobject  is 
composed  of  polygons  or  it  is  a  sphere.  Using  "igure  2.3a  as  an  example  again, 
the  left  and  right  weights  are  spheres  and  the  bar,  instead  of  being  a  perfect 
cylinder,  is  composed  of  polygons  and  actually  has  an  octagonal  shape,  Figure 
2.3b.  Aside  from  containing  a  pointer  to  the  common  part  record,  examined  next, 
and  data  for  its  bounding  volume,  it  also  contains  information  on  the  subobject 
type.  This  subobject  type  field  indicates  the  geometry  of  the  subobject,  i.e.,  it  is 
either  a  sphere  or  a  polygonal  object,  which  is  an  object  composed  of  polygons. 
?h>  r.:or::in:;on  >  't^red  .e<-a:is«*  urferenr  r.  routines  .ire  useu  tor  -acn. 

object  type.  Currently  a  1  indicates  that  planar  intersection  routines  should  be 
used  and  a  0  indicates  spherical  intersection  routines  should  be  used. 


A 


5.  Common  Parts 


A  good  illustration  of  common  parts  is  found  in  an  ordinary 
checkerboard.  Figure  2.4a.  In  this  figure  one  object  exists— the  checkerboard.  This 
in  turn  has  one  subobject,  itself.  This  subobject  has  two  common  parts:  the  white 
squares  and  the  black  squares.  Table  2.1  listed  fhe  characteristics  of  an  object 
and  it  is  in  the  common  parts  record  that  these  characteristics  are  stored.  Each  of 
these  common  parts  records  contains  a  pointer  to  an  array  of  polygon  records.  It 
is  through  this  arrangement  that  the  polygons  inherit  the  characteristics  of  the 
object.  Therefore,  the  common  parts  array,  also  called  the  Cparts  array,  is  an 
array  of  records,  with  one  array  per  subobject,  and  each  common  parts  record 
points  to  its  own  set  of  polygons,  Figure  2.4b. 

6.  Polygons 

The  polygons  array,  too,  is  an  array  of  records  with  one  array  for  each 
subobject.  This  is  the  smallest  physical  item  in  the  scene  and  the  one  against 
which  the  actual  intersections  are  determined. 

7.  Vertex  Array 

The  vertex  array  is  an  array  of  points  that  define  the  polygons  that 
compose  the  subobject. 

This  data  structure  as  presented  by  Falby  [Ref.  3]  proved  itself  to  be  both 
flexible  and  easy  to  use.  An  example  of  a  data  base  that  used  this  structure  and 
which  was  used  in  testing  this  ray  tracer  can  be  seen  in  Appendix  B. 


Figure  2.4a  -  Example  of  Object  with  Two  Common  Parts 


Figure  2.4b  -  Record  Layout  of  2.4a 


III.  RAY  TRACING  INTERSECTION  CONSIDERATIONS 


The  methodology  behind  ray  tracing  is  quite  simple.  However,  it  does  require 
an  understanding  of  the  fundamentals  of  vector  calculus  and  geometric  optics.  A 
discussion  of  the  fundamentals  is  beyond  the  scope  of  this  study.  An  overview  of 
some  of  the  fundamentals  is  in  order. 

A.  RAY  TRACING  MECHANICS 

By  its  very  definition,  ray  tracing  is  simply  the  tracing,  or  following,  of  a  ray 
from  its  source  through  space  and  determining  any  possible  intersections  that  may 
occur  between  it  and  an  object.  The  natural  way  to  model  a  ray  in  order  to  do 
this  tracing  is  by  using  vectors.  A  vector  is  not  only  a  precise  way  to  represent  a 
ray  but  the  basic  operations  on  vectors  in  three  space,  addition,  subtraction,  dot 
product,  and  cross  product  provide  the  tools  necessary  to  determine  the 
intersections.  These  tools,  along  with  other  techniques  found  in  vector  calculus 
and  geometric,  optics  provide  the  means  to  deal  with  the  two  problems 
encountered  in  ray  tracing,  i.e.,  the  ray  direction  determination  problem  ans  the 
intersection  problem. 


1.  The  Ray  Direction  Problem 


position,  usually  referred  to  as  the  view  ray,  is  the  simplest  to  solve.  Every  point 
in  a  coordinate  system  can  be  associated  with  a  ray,  and  determining  the  direction 


of  a  ray  between  two  points  can  be  solved  by  using  vector  subtraction.  The  last 
step  in  the  ray  tracing  process  is  determining  what  takes  place  when  a  ray 
intersects  an  object.  This  requires  the  application  of  the  laws  of  geometric  optics. 
Once  a  ray  strikes  an  object,  either  one  or  two  additional  rays  will  be  generated. 
These  new  rays  are  referred  to  as  the  reflected  and  refracted  rays,  Figure  3.1.  The 
three  basic  laws  of  reflection  and  refraction  are  listed  as  [Ref.  8:  pp:  32-33] 

1.  The  incident,  reflected,  and  transmitted  rays  all  reside  in  a  plane,  known 
as  the  plane  of  incidence,  which  is  normal  to  the  surface  of  the  object. 

2.  The  angle  of  incidence  is  equal  to  the  angle  of  reflection  ©i  =  ©r. 

3.  The  incident  and  transmitted  ray  directions  are  related  by  Snells’  law: 
nsin@t  =  ntsin©(. 


An  illustration  of  these  laws  is  shown  in  Figure  3.1.  Rogers  [Ref.  1:  p.  367] 
provides  a  method  for  determining  the  direction  of  the  reflected  and  refracted 
rays.  The  direction  of  r,  the  reflection  ray,  and  p,  the  refraction  ray  are  given  as: 


where 


r  =  v'  +2  n 


n 


v 

t 

v  =  - 

vri 
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kf  =  (*„V|2  ~  K  +  n\2)  2 

r?2 

kn  =  ~ 

'h 

Here  A:  is  the  ratio  of  refractive  indices,  k,  is  the  Fresnel  coefficient,  v  is  the 
normal  vector  in  the  direction  of  the  incoming  ray,  v'  is  the  unit  normal  vector  in 
the  direction  of  the  incoming  ray,  n  is  the  unit  surface  normal,  and  rj}  and  rj2  are 
the  refraction  coefficients  for  mediums  the  rays  pass  through.  This  is  illustrated 
in  Figure  3.2. 

2.  The  Intersection  Problem 

As  stated  above  the  intersection  computation  is  the  most  time 
consuming  part  of  the  ray  tracing  process.  It  is  not  that  the  process  itself  is  so 
difficult  but  because  several  steps  need  to  be  done  for  each  iteration.  Two  types 
of  intersection  computations  are  required  to  be  performed:  determining  the 
intersection  between  a  line  and  a  sphere  and  determining  the  intersection  between 
a  line  and  a  polygon.  The  first  type  is  the  simplest  to  solve  and  is  why  the  sphere 
is  generally  used  as  the  bounding  volume.  The  calculation  of  the  intersection 
point  between  a  line  and  a  sphere  involves  solving  the  equation  for  the  line  and 
■  r.e  ->oner»*  -lniuitaneousiy.  The  sphere  -s  aerined  >v  me  -quurion 

(j:  -  a)2  +  (y  -  j3)2  +  (z  -  6)2  =  2  (3.1) 

where  (a.  /?,  S)  is  the  center  point,  r  is  the  radius,  and  (r,  y,  z )  is  a  point  on  the 
sphere.  The  line  is  defined  by  the  parametric  equations 
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x  =  at  +  xQ  y  =  bt  +  yQ  z  =  ct  +  zQ  (3.2) 

where  ( xQ ,  y0,  20)  is  a  known  point  on  the  line,  and  a,  b,  and  c  are  coefficients 
from  *he  vector  of  -  hj  —  ck  which  is  parallel  to  the  line.  These  equations  must 
first  be  solved  for  t.  The  solution  to  this  provides  two  answers.  First,  it  indicates 
whether  or  not  an  intersection  actually  takes  place.  Second,  if  there  is  one,  it 
indicates  how  many  intersections,  either  one,  in  the  case  where  the  line  is  tangent 
to  the  sphere,  or  two,  in  the  case  where  it  actually  enters  the  sphere.  In  the  case 
of  two  intersection  points,  a  check  must  be  done  to  determine  which  is  closer  to 
the  origin  of  the  ray. 

The  intersection  between  a  line  and  a  polygon  is  more  involved.  This 
problem  is  comprised  of  two  parts:  determining  the  intersection  point  between  a 
line  and  a  plane,  and  determining  whether  or  not  the  intersection  point  lies  within 
the  polygon.  Like  the  line-sphere  intersection  problem,  this  one  also  involves  the 
solving  of  two  equations  simultaneously.  The  first  of  these  is  the  equation  for  a 
plane  which  is  defined  as 

Ax  +  By  +  Cz  =  D  (3.3) 

where  A,  B,  C,  and  D  are  constants  and  (x,  y,  z)  is  a  point  on  the  plane.  The 

second  equation  used  is  'he  parametric  equation  that  defines  a  iine.  So.  3.2.  The 

solution  of  this  requires  first  substituting  the  equations  for  the  line  into  the 

equation  of  the  plane.  The  result  is  an  equation  in  t,  which  when  solved  and 

substituted  back  into  the  equations  for  the  line,  provides  the  intersection  point 
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between  the  line  and  the  plane.  Once  this  point  is  calculated,  it  is  then  necessary 
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to  determine  whether  or  not  it  lies  within  the  polygon.  A  general  way  of  solving 
this  is  to  simply  determine  the  relationship  between  the  intersection  point  and 
each  edge  of  the  polygon.  The  point  that  lies  on  the  inside  of  each  edge  also  lies 
within  the  polygon.  If  it  fails  the  test  for  any  edge,  then  it  lies  outside  the 
polygon.  The  drawback  to  this  approach,  is  that  it  only  works  for  convex 
polygons.  In  this  study  we  assume  all  polygons  are  convex.  [Ref.  9], 

B.  THE  RAY  DATA  STRUCTURE 

Rogers  [Ref.  l]  suggests  a  data  structure  for  a  ray  in  a  ray  tracer.  It  is  that 
data  structure  which  is  used  in  this  study.  Table  3.1  lists  the  data  used  to  model 
each  ray.  We  examine  each  item  of  this  structure  as  adapted  from  Rogers  [Ref.  1: 
p.  373]. 


FIELD  NAME 
Ray  Type 
Ray  Origin 
Ray  Vector 
Source  Ray  Type 
Intersection  Flag 
Object  Index 
Subobject  Index 


TABLE  3.1  -  RAY  DATA 

VARIABLE  NAME 
type 
*»  !/>  ^ 

*»  y, 2 

Stype 

flag 

obj  idx 

subobj  :dx 


VALUE 

enumerated  or  coded 
real 
real 

enumerated  or  coded 
boolean  or  coded 
integer 
:nteger 


Common  Part  index 

'  .'part  idx 

nroger 

Polygon  index 

polygon  idx 

integer 

Intersection  Point 

y,  2 

real 

Distance 

d 

real 

Transmitted  Intensity 

It 

real 

Specular  Intensity 

Is 

real 

1.  Ray  Type 


< 

i 


The  ray  type  field  identifies  a  ray  as  either  a  view  ray,  a  reflected  ray,  or 
a  refracted  ray.  The  values  put  in  this  field  are  generally  of  an  enumerated  type 
and  consist  of  reflected,  refracted,  view,  and  none. 

2.  Ray  Origin 

The  ray  origin  field  contains  the  point  that  identifies  the  position  from 
which  the  ray  originated.  For  instance,  if  it  is  the  view  ray,  its  point  of  origin  is 
the  view  position.  If  it  is  a  reflected  or  refracted  ray,  its  origin  is  the  intersection 
point  that  it  originated  from. 

3.  Ray  Vector 

The  ray  vector  field  contains  the  vector  heading  of  the  ray. 

4.  Source  Ray  Type 

The  source  ray  type  field  contains  the  ray  type  of  the  source  ray  for  this 
particular  ray.  For  instance,  a  view  ray  does  not  have  any  source  ray  as  it  is  the 
starting  ray  for  the  process.  Hence,  none  is  in  the  type  field.  If  the  view  ray 
intersects  an  object  and  both  a  reflected  and  refracted  ray  are  generated,  then  the 
source  ray  for  both  of  them  is  the  view  ray.  Likewise  if  the  reflected  or  refracted 
ray  hits  an  object  and  creates  further  rays  then  it  becomes  the  source  ray  for 
rhose  rays  in  creates. 

5.  Intersection  Flag 

Originally  the  intersection  flag  is  set  to  false  and  it  is  only  set  to  true 

when  there  is  an  intersection  between  this  ray  and  an  object. 
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6.  Object  Index 


The  object  index  provides  an  index  into  the  array  of  object  records 
making  it  possible  to  select  any  object  easily. 

7.  Subobject  Index 

The  subobject  index  provides  an  index  into  the  array  of  subobject 
records  and  helps  uniquely  identify  each  subobject. 

8.  Common  Part  Index 

The  common  part  index  provides  the  index  into  the  array  of  common 
part  records  uniquely  identifying  each  common  part  record. 

9.  Polygon  Index 

The  polygon  index  provides  the  index  into  the  array  of  polygon  records 
and  uniquely  identifies  each  polygon. 

10.  Intersection  Point 

The  intersection  point  field  holds  the  position  of  the  intersection  point 
between  the  current  ray  and  an  object. 

11.  Distance 

The  distance  field  contains  the  distance  between  the  current  ray’s  point 
of  origin  and  its  point  of  intersection. 


Transmitted  Intensity 


The  transmitted  intensity  field  contains  the  red,  green,  and  blue 
intensity  values,  in  a  range  between  0  and  1,  of  the  light  that  is  incoming  along 

any  refracted  ray  that  this  ray  produces. 
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13.  Specular  Intensity 

The  specular  intensity  field  contains  the  red,  green,  and  blue  intensity 
values,  in  a  range  between  0  and  1,  of  the  light  that  is  incoming  along  any 
specular  ray  that  this  ray  produces. 

C.  INTERSECTION  METHODOLOGY 
1.  Intersecting  a  Planar  Polygon 

a.  Generating  the  Initial  Ray 

The  generation  of  the  initial  or  view  ray  is  shown  in  Figure  3.3. 
This  ray,  p  -  v,  is  calculated  using  vector  subtraction.  The  two  vectors  used  are 
the  ones  associated  with  the  points  for  the  view  position,  v,  and  the  pixel  through 
which  the  ray  passes,  p.  The  ray  associated  with  the  view  position  is  to  be 
subtracted  from  the  ray  associated  with  the  pixel. 

b.  Intersecting  the  Bounding  Volumes 

After  each  ray  is  generated,  each  object  in  the  object  list  is  checked, 
one  at  a  time,  to  determine  whether  or  not  the  ray  strikes  any  of  the  bounding 
containers  of  the  objects  in  the  scene.  In  our  implementation,  the  bounding 
container  is  the  sphere,  which  is  reduced  to  a  bounding  circle,  C,  see  Figure  3.4. 
This  ••jrcie‘»  radius  .s  "lie  -ame  as  me  raaius  of  me  bounding  sphere.  .V .  and  :t  .ies 
on  the  plane,  P ,  which  contains  the  center  point,  q,  of  the  bounding  sphere.  The 
inverse  of  the  incoming  ray,  tJ,  is  the  surface  normal  of  this  plane.  It  is  not 
necessary  to  determine  where  on  the  bounding  sphere  a  ray  hits  since  at  this  stage 


we  are  only  interested  in  just  finding  out  whether  or  not  it  hits  it.  Because  we  are 


not  concerned  with  where  the  ray  strikes  the  bounding  sphere,  a  bounding  circle  is 
used.  Determining  the  intersection  with  a  circle  requires  less  work  than 
determining  the  intersection  with  a  sphere. 

The  first  step  in  tracing  the  ray  begins  by  taking  each  object  and 
constructing  a  bounding  circle.  Once  this  circle  is  constructed,  the  intersection 
point  of  the  view  ray  with  it  is  calculated.  Then  the  distance  between  the 
intersection  point  and  the  center  point  is  determined  and  compared  to  the  radius 
of  the  bounding  circle  to  see  if  it  falls  within  the  circle.  If  it  does  fall  within  the 
circle,  this  indicates  that  the  ray  intersected  the  bounding  volume.  This  process 
then  needs  to  be  repeated  for  each  object.  If  an  object’s  bounding  circle  is  hit, 
this  process  must  then  be  repeated  for  each  subobject  of  that  object, 
c.  Intersecting  the  Polygon 

Once  a  particular  subobject  is  identified  by  the  bounding  volume 
tests,  the  common  parts  list  is  processed,  one  record  at  a  time.  Each  one  of  these 
common  parts  records  contains  a  pointer  to  those  polygons  that  make  up  the 
object  being  rendered.  This  list  of  polygon  records  is  then  processed  after  the 
bounding  volume  processing  is  completed.  The  entire  list  needs  only  to  be 
processed  anui  an  intersection  :s  round.  The  processing  involved  in  rhis  s  'he 
most  computationally  expensive  part  of  the  entire  ray  tracing  algorithm.  This 
computation  consists  of  three  steps:  determining  the  orientation  of  the  polygon, 
calculating  the  intersection  point  between  the  ray  and  the  plane  that  contains  the 


polygon,  and  determining  whether  or  not  the  intersection  point  lies  on  the 
polygon. 

Ill  Establishing  the  Orientation.  First,  each  polygon  needs  to  be 
checked  for  a  correct  orientation.  This  is  a  straightforward  step  carried  out  by 
calculating  the  angle  between  the  surface  normal  of  that  polygon  and  the  inverse 
of  the  view  ray.  If  it  is  90  degrees  or  greater,  it  is  facing  in  the  wrong  direction  to 
be  intersected.  If  it  is  less  than  90  degrees,  the  next  step  is  to  determine  whether 
the  ray  intersects  the  particular  plane  that  that  particular  polygon  lies  on.  The 
first  step  in  doing  this  is  to  determine  the  equation  for  the  plane  in  which  that 
polygon  lies. 

(2)  Intersecting  a  Plane.  If  the  correct  orientation  exists  for  a 
polygon  to  be  intersected,  the  next  step  is  to  define  the  plane  containing  the 
polygon  of  interest.  The  equation  for  a  plane  was  given  earlier  as 

Ax  -f-  By  +  Cz  =  D 

where  A,  B,  C,  and  D  are  constants  and  can  be  calculated  by  the  following 
equation. 

-  Vl(Z2  ~  Z3^  ~  V2{ZS  ~  zi)  *  iM*!  ~  z2) 

3  -  :. 

c  =  *i(v2  -  ^  -  y.)  *  *1(1/1  -  v2) 

d  =  ~z\iy2zz  ~  v2z2)  -  x2  (1/3^1  *  V\ZJ  ~  1 3  (1/1*2  i/»~i) 

where  (jtj,  yv  2J,  (xv  yv  z2),  and  (xy  yv  z^)  are  points  of  the  vertices  of  the 


polygon.  Once  the  equation  for  the  plane  is  known,  it  must  then  be  solved 
simultaneously  with  the  equation  of  the  line  representing  the  ray. 


(3)  Location  of  a  point  with  respect  to  a  polygon.  Once  the 


intersection  point  is  determined,  it  needs  to  be  checked  to  see  whether  or  not  it 
lies  within  the  polygon.  This  is  the  most  computationally  expensive  part  of  the 
process.  This  process,  see  Figure  3.5,  requires  that  a  plane,  P,  called  the  bounding 
plane,  containing  two  vertices,  for  example  A  and  B,  of  an  edge  and  another 
arbitrary  point,  A',  not  on  the  polygon  or  its  plane,  be  constructed.  This  must  be 
done  for  each  edge.  Once  the  bounding  plane  is  constructed,  the  point  in  question 
must  then  be  checked  to  see  whether  or  not  it  lies  to  the  polygon’s  side  of  the 
plane  that  now  contains  the  edge  of  the  polygon.  This  is  done  by  plugging  in 
another  vertex  of  the  polygon  into  the  equation  for  the  plane  that  was  just 
constructed,  and  then  plugging  in  the  intersection  point.  If  the  results  from  these 
two  equations  have  the  same  relationship,  i.e.  if  the  sign  is  the  same,  then  the 
intersection  point  lies  on  the  polygon  side  of  the  bounding  plane.  This  check 
must  then  be  repeated  for  each  edge  of  the  polygon.  In  order  for  the  intersection 
point  to  lie  on  the  face  itself,  it  must  be  found  to  lie  on  the  polygon’s  side  of  each 
»dge  .V  ^ h is  ooint.  it  is  determined  whether  or  not  the  ray  strikes  the  object.  If 
\.e  '  :iv  >nje<-r.  -he  pointers  *o  'ius  polygon  must,  be  stored  .n  '.he 

ray  structure  along  with  the  intersection  point.  This  information  is  needed  later 
on  when  determining  the  direction  of  reflected  and  refracted  rays  as  well  as  the 


intensity  of  the  light  that  is  reflected  from  this  position. 
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2.  Intersection  of  a  Sphere 


The  sphere  is  the  easiest  object  to  work  with  in  a  ray  tracer.  Since  the 
sphere  can  act  as  its  own  bounding  volume,  the  center  point  and  radius  are 
already  available  in  the  subobject  record.  This  eliminates  the  need  for  the  polygon 
array.  The  center  point  and  radius  are  the  only  information  necessary  to  model  a 
sphere  for  a  ray  tracer.  Determining  the  intersection  of  a  line  and  a  sphere  is 
nothing  more  than  the  simultaneous  solving  of  their  equations  for  the  variable  t  . 
The  solution  to  this  gives  a  quadratic  equation  in  t  which  can  then  easily  be 


solved. 


IV.  THE  INTENSITY  PROBLEM 


One  of  the  strong  points  of  the  ray  tracing  algorithm  is  that  a  global 
illumination  model  can  easily  be  integrated  into  it.  In  fact,  ray  tracing  and  global 
illumintion  models  seem  to  naturally  complement  each  other.  A  global 
illumination  model  takes  into  account  all  of  the  light  sources  in  the  scene  in 
calculating  the  intensity  at  each  point.  This  means  taking  into  account  ♦he 
ambient  light  that  exists  in  a  scene,  light  that  comes  directly  from  a  particular 
light  source(s),  and  light  that  is  reflected  off  an  object(s).  It  also  includes  the 
coefficients  necessary  to  model  the  way  an  object  reacts  with  light.  A  great  deal 
of  work  has  been  done  in  this  area.  The  most  notable  model  is  the  VVhitted 
illumination  model,  and  it  is  the  one  that  has  been  implemented  here  [Ref.  1:  pp. 
365-366].  The  Whitted  algorithm  is  based  on  the  three  models  shown  in  Figures 
4.1,  4.2,  and  4.3.  These  models  will  now  be  examined  more  closely.* 

A.  LOCAL  ILLUMINATION  MODELS 
1.  Diffuse  Reflection  Model 

The  ffrst  of  rhese  models  is  a  perfect  diffuser.  Such  i  model  :s  governed 

by  Lambert’s  cosine  law.  This  law  states  that  the  intensity  of  light  reflected  from 
a  perfect  diffuser  is  proportional  to  the  cosine  of  the  angle  between  the  light 
‘The  contents  of  this  chapter  are  close  adaptations  from  Rogers  [Ref.  ij  and  Falby  [Ref.  3j. 


direction  and  the  normal  to  the  surface.  This  can  be  expressed  mathematically  as 
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where  1  ;s  the  reriected  intensity,  *s  the  incident  .ntensity  from  a  point  .ight 
source,  kd  is  the  diffuse  reflection  constant,  unique  to  each  object,  and  ©  is  the 
angle  between  the  light  direction  and  the  surface  normal,  see  Figure  4.1.  Since  the 


diffuse  reflection  coefficient  kd  varies  from  material  to  material  and  is  also  a 


function  of  the  wavelength  of  the  light,  it  is  often  easier  to  just  assume  it  a 

constant  for  simole  illumination  models.  rRef.  2:  o.  3121 


2.  Specular  Reflection  Model 


The  second  model  illustrates  the  characteristics  of  specular  reflection 
which  is  directional,  unlike  diffuse  reflection.  This  means  that  the  greatest 
intensity  of  the  specularly  reflected  light  can  only  be  seen  if  the  view  angle 
coincides  with  the  reflection  angle,  Figure  4.2.  The  further  off  the  viewing  angle  is 
from  the  reflection  angle,  the  dimmer  the  intensity  becomes.  Because  of  the 


complex  physical  characteristics  of  specularly  reflected  light,  an  empirical  model 
due  to  Bui-Tuong  Phong  is  usually  used  for  simple  illumination  models  [Ref.  10]. 
This  is  expressed  mathematically  as 


X’l  (-)  |  cos  '< 


where  w(&  ,  A),  the  reflection  curve,  gives  the  ratio  of  the  specularly  reflected 
light  to  the  incident  light  as  a  function  of  the  incidence  angle  0(  and  the 


wavelength  A.  Because  w(0(,  A)  is  such  a  complex  function,  it  is  frequently 
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replaced  by  an  aesthetically  or  experimentally  determined  constant  k}  which  then 
yields 

I  —  r,k  cosna 

S  4  J 

Also,  n  is  a  power  that  approximates  the  spatial  distribution  of  the  specularly 
reflected  light.  Typically  a  value  of  200  for  n  is  used  to  model  a  very  shiny 
surface  and  a  value  of  10  is  used  for  a  dull  surface  [Ref.  3:  p.  72].  [Ref.  1:  pp.  313- 
315] 

3.  Combined  Model 

If  just  point  sources  are  assumed,  as  in  the  two  models  just  discussed, 
any  object  not  receiving  light  directly  from  the  source  appears  black.  In  order  to 
properly  render  a  scene,  it  is  also  necessary  to  take  into  account  ambient  light,  the 
light  that  is  reflected  off  other  surfaces.  Including  a  model  for  ambient  light  into 
the  intensity  calculations  is  not  feasible.  Ambient  light  represents  a  distributed 
light  source  and  as  such  is  a  very  complex  function.  Therefore  it  is  treated  as  a 
constant  diffuse  term  and  linearly  combined  with  the  other  terms.  Also  not 
included  in  the  above  model  is  the  effect  that  distance  has  on  light.  It  is  well 
known  that  the  farther  away  an  object  or  light  source  is,  the  dimmer  it  gets.  The 

actual  tormina  'o  produce  that  proper  attenuation  is  — .  where  rhe  intensity  of 

of 

light  decreases  as  the  square  of  the  distance  from  the  source  increases.  However  it 
has  been  shown  that  linear  attenuation  can  actually  produce  more  realistic  results. 
With  these  two  additions  the  complete  model  now  looks  like: 
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where  Ia  is  the  incident  ambient  light  intensity,  ka  is  the  ambient  diffuse 
reflection  constant  and  K  is  an  arbitrary  constant  that  can  be  used  to  adjust  the 
results.  [Ref.  1:  p.  313] 

The  complete  model  just  presented  can  now  be  modified  to  better  fit  in 
with  a  ray  tracer.  Recalling  the  formula  for  the  dot  product  of  two  vectors  allows 
writing  the  cos©  as  n-L  and  writing  cosa  a s  R-S  which  gives  us: 

h 

I  =  hka  +  - [^(rf-I)  +  *,(jf-sn 

d  +  K 


So  far  we  have  only  examined  the  case  where  just  one  light  source  is  present.  If 
there  are  several  light  sources,  the  effects  are  added  linearly,  and  the  equation  now 
becomes: 


i  =  KK  +  S — —  [*„(*•£,-)  +  *.(*,■*)“] 

J  =  1  + 

This  then  is  the  complete  local  illumination  model.  [Ref.  1:  pp.  312-316] 


B.  GLOBAL  ILLUMINATION  MODEL 

The  ••.ompiete  local  Humiliation  modei  just  presented  forms  the  oasis  for  the 

algorithm  that  was  implemented  for  this  study  [Ref.  1:  pp.  363-378]  see  Figure 
4.3. 


'  =  Etj  -i '  £,]  +  *,Et,(  S  R, I  +  K1,  +  V, 

/  ; 

In  the  above  equation  &a,  kd,  ks ,  and  kt  are  the  ambient,  diffuse,  specular 
reflection,  and  transmission  coefficients,  all  of  which  have  values  between  0  and  1. 
Jo,  Is,  /(,  and  7;  are  the  intensities  of  the  ambient  light,  the  specularly  reflected 

light,  the  transmitted  light,  and  the  light  directly  from  a  light  source.  These  also 
hold  values  between  0  and  1.  The  remaining  variables  n,  L .,  S,  and  are  the 
surface  normal  at  the  intersection  point,  the  direction  of  the  jth  light  source,  the 
local  sight  vector,  and  the  local  reflection  vector  from  the  jth  light  source.  A 
careful  comparison  between  this  model  and  the  complete  local  illumination  model 
reveals  that  the  only  new  terms  are  the  J^and  It  terms.  These  are  the  terms  used 
to  account  for  the  light  that  comes  in  along  the  reflected  and  refracted  rays  that 
originated  at  this  point.  If  holds  the  intensity  for  the  reflected  ray  and  If  holds 
the  intensity  for  the  refracted  ray.  These  two  values  in  turn  are  calculated  using 
the  exact  same  model.  For  the  last  intersection  point  in  the  scene,  the  one  whose 
reflected  and  refracted  rays  do  not  intersect  anything,  Ia  and  It  are  set  to  0.  The 
k  and  k(  terms  are  coefficients  included  to  better  model  the  way  this  object  reacts 
with  the  light  incoming  along  the  reflected  and  refracted  rays. 

This  'hen  is  the  compiet*’  giobai  illumination  model  used  in  this  study.  It  is 
the  simplicity  of  this  algorithm  that  makes  it  so  easy  to  understand  and 
implement.  In  essence,  it  is  saying  that  the  output  intensity  is  nothing  more  than 
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a  sum  of  all  possible  light  sources  with  the  coefficients  determining  the  intensity  of 
light  that  comes  from  a  particular  object. 


►  v 


■  / 


V.  RAY  TRACING  ALGORITHM 


Rogers  suggested  an  algorithm  for  ray  tracing  [Ref.  1:  pp.  374-377).  It  is  that 


algorithm  that  has  been  the  basis  for  this  study.  The  following  is  a  description  of 


that  algorithm  as  it  has  been  implemented  here. 


A.  TRACING  THE  RAYS 


To  begin  the  ray  tracing  process,  the  first  thing  done  is  the  determination  of 


the  direction  of  the  view  ray.*  The  ray  data,  mentioned  in  chapter  III,  is  then 


initialized.  After  the  view  ray  is  generated  and  the  ray  data  initialized,  the  ray, 


which  is  represented  by  this  ray  data,  is  pushed  onto  the  stack,  which  is  used  to 


implement  the  ray  tracing  tree.  The  process  then  moves  into  the  actual  ray 


tracing  loop. 


Once  in  the  loop,  the  stack  is  checked  to  see  if  it  is  empty.  If  it  is  not  empty, 


the  stack  is  then  popped  to  access  the  ray  information.  The  first  thing  checked  is 


whether  or  not  the  intersection  flag  is  set.  If  it  is,  that  indicates  that  that 


particular  ray  has  already  been  terminated  (by  hitting  an  object),  and  the  process 


if  determining  "he  intensities  begins. 


‘Each  ray,  is  modeled  as  a  vector,  and  is  converted  to  a  unit  vector  immediately  after  its 
determination. 


If  the  intersection  flag  is  not  set,  a  new  ray  has  been  generated,  either  a  new 
view  ray  or  a  new  reflected  or  refracted  ray.  All  of  these  rays  are  grouped  under 
the  more  general  title  of  a  shooting  ray.  At  this  stage,  the  rav  must  be  sent 
through  the  intersection  procedures  to  determine  w'hether  or  not  it  intersects  any 
object.  The  intersection  routines  start  at  the  highest  level  of  the  picture  and  step 
through  the  linked  list  of  objects,  subobjects,  and  common  parts  to  the  actual 
polygons. 

At  the  top  level  of  the  intersection  routines,  each  one  of  the  objects  is 
checked.  First,  they  are  checked  to  see  if  there  is  an  intersection  with  the  object’s 
bounding  circle.  If  there  is  an  intersection,  the  distance  between  the  intersection 
point  and  the  origin  of  the  ray  is  calculated  and  placed  in  the  ray  data.  Second,  a 
check  is  done  to  insure  that  the  objects  lie  in  front  of  the  ray’s  origin.  In  front 
refers  to  those  objects  that  lie  in  the  direction  of  the  shooting  ray  from  its  point  of 
origin.  Since  the  ray  is  being  modeled  by  a  vector,  which  only  indicates  direction, 
every  object  along  the  line  described  by  the  vector  and  the  origin  of  the  ray  is 
considered,  see  Figure  5.1.  The  way  to  test  for  this  is  to  generate  a  test  vector 
between  the  origin  of  the  ray  and  any  intersection  points  of  the  shooting  ray, 
eliminating  the  intersection  ooints  that  lie  in  the  oooosite  direction  from  the 
mooring  .‘ay. 

As  each  object  is  processed,  the  objects  whose  bounding  volumes  are 
intersected  by  the  shooting  ray  are  processed  further  to  see  if  any  of  their 
subobjects  are  also  intersected.  The  same  basic  procedure  for  finding  an 


intersection  is  used  here  as  is  used  for  finding  an  intersection  for  the  objects.  If  the 
object  and  the  subobject  are  one  and  the  same,  this  check  need  not  be  done. 

For  each  subobject  whose  bounding  volume  is  intersected,  the  next  step  is  to 
check  the  type  of  geometric  object,  i.e.,  a  sphere  or  a  polygonal  object.  This  needs 
to  be  done  in  order  to  determine  which  intersection  routines  are  needed.  In  the 
case  of  a  sphere,  it  becomes  just  a  matter  of  solving  two  simultaneous  equations, 
discussed  earlier— one  for  a  line,  the  other  for  the  sphere.  When  dealing  with  a 
polygonal  object,  it  becomes  more  complicated  because  each  face  of  the  polygon 
must  be  checked.  The  first  thing  that  needs  to  be  done  is  to  check  the  orientation 
of  the  face  in  question.  When  the  polygon  does  have  the  right  orientation,  the 
intersection  process  continues.  The  first  step  is  to  determine  the  equation  of  the 
plane  that  contains  the  polygon.  This  plane  is  calculated  from  any  three  vertices 
of  the  face  in  question.  Once  the  plane  has  been  established,  the  intersection  point 
between  it  and  the  shooting  ray  is  calculated.  With  the  intersection  point  thus 
established,  the  next  step  is  to  determine  whether  or  not  the  point  lies  on  the  face 
of  the  polygon.  If  the  intersection  point  is  found  to  lie  on  the  face  of  the  polygon, 
then  this  point  is  saved  and  placed  in  the  ray  data. 

If  no  intersection  is  found  and  the  shooting  ray  is  either  a  reflecting  or 
r-'tr*ic':ng  ray.  ' r  >  uisc-imeu  mu  :;»>  n***!>irv  a.c u.a r :< >n  -if ><•■'-  » .: 
the  shooting  ray  is  a  view  ray,  then  the  intensity  is  set  to  the  background 
intensity  which  is  then  displayed  and  the  next  ray  is  generated 
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When  there  is  an  intersection,  the  first  thing  to  be  checked  is  whether  or  not 
there  is  enough  room  on  the  stack.  Since  the  stack  holds  only  part  of  the  ray  tree 
at  any  one  time,  it  need  only  be  long  enough  to  contain  the  longest  anticipated 
branch.  A  particular  branch  of  the  ray  tree  is  terminated  when  both  the  reflected 
and  refracted  rays  at  an  object  intersection  leave  the  scene  or  when  the  available 
stack  length  is  exceeded.  When  both  rays  leave  the  scene,  their  contribution  to 
the  illumination  at  the  source  ray  is  zero.  W’hen  the  available  stack  length  is 
exceeded,  the  algorithm  calculates  the  illumination  at  the  source  ray  using  only 
the  ambient,  diffuse,  and  specular  reflection  components  at  the  source  ray 
intersection.  An  extension  algorithm  is  given  in  Rogers  whereby  the  algorithm  can 
be  extended  one  additional  depth  in  the  tree  without  exceeding  the  maximum 
stack  depth.  However,  the  implementation  of  that  was  not  necessary  here.  When 
the  stack  does  get  full,  it  becomes  a  matter  of  calculating  the  intensity  at  that 
intersection  point  and  setting  the  appropriate  value  I(  or  in  the  source  ray. 
Ref.  1 :  p.  372] 

When  the  stack  is  not  full,  then  the  distance  between  the  source  point  of  the 
shooting  ray  and  the  intersection  point  is  determined  and  placed  in  the  ray  data. 
The  rav  is  then  placed  hark  on  the  stack.  Once  that  is  done,  any  reflecting  and/or 
■■  "r  !<  .1-4  "iv>  ri.ana  : 4  ;:«»  ntersee-iou  < n:it  le'emune-i.  "tv 

data  initialized,  and  then  placed  on  the  stack  as  new  rays  with  the  reflecting  ray 
being  placed  on  first.  It  is  important  to  keep  this  order  of  rays  in  rnind  because  it 
is  necessary  to  know  the  number  of  rays  to  pop  when  setting  the  intensities  of  the 


source  rays.  With  these  new  rays  in  place  on  the  stack,  the  program  returns  to 


the  beginning  of  the  ray  tracing  cycle.  In  the  absence  of  reflecting  or  refracting 
rays,  then  the  first  rav  popped  is  the  view  ray.  Now  since  this  ray  already  has  its 
intersection  flag  set,  the  intensity  at  the  point  of  its  intersection  is  calculated  and 
displayed.  If  the  ray  popped  is  a  reflective  or  refractive  ray,  then  it  becomes  the 
new  shooting  ray.  This  cycle  continues  until  either  no  more  reflecting  or 
refracting  rays  are  produced  or  until  the  stack  becomes  full.  This  process  is 
summarized  by  the  pseudocode  description  in  Figure  5.2. 

B.  DETERMINING  THE  INTENSITY 

If  a  ray’s  intersection  flag  is  found  set  at  the  beginning  of  the  ray  tracing 
process,  that  ray  is  sent  into  the  intensity  processing  part  of  the  ray  tracer.  The 
first  step  in  this  process  is  to  take  the  ray  data  and  to  determine  the  intensity  at 
that  ray’s  intersection  point.  This  intensity  is  then  divided  by  the  distance 
between  the  ray’s  point  of  origin  and  its  intersection  point  in  order  to  properly 
attenuate  it.  This  process  is  done  for  each  of  the  primary  colors  —  red,  green,  and 
blue.  If  a  view  ray  is  being  considered,  that  means  that  it  was  the  last  ray  on  the 
stack.  Therefore,  it  is  the  final  intensity  to  be  calculated  and  hence  is  the  actual 


READ  IN  DATA  FILE 


for  Y  1  to  MAX  ROWS  io 

for  X  :=  1  to  MAXCOLUMNS  do 
INITIALIZE  VIEW  RAY 
PUSH  RAY  ONTO  STACK 
repeat 

POP  RAY  FROM  STACK 
if  INTERSECTION  FLAG  SET  then 
CALCULATE  INTENSITY 
else 

CHECK  FOR  INTERSECTION 
if  INTERSECTION  FLAG  SET  then 
if  STACK  EXCEEDED  then 
CALCULATE  INTENSITY 
else 

PUSH  RAY  BACK  ON  STACK 
CALCULATE  REFLECTED  RAY 
CALCULATE  REFRACTED  RAY 
if  REFLECTED  RAY  EXISTS  then 
INITIALIZE  REFLECTED  RAY 
PUSH  REFLECTED  RAY  ON  STACK 
end  if 

if  REFRACTED  RAY  EXISTS  then 
INITIALIZE  REFRACTED  RAY 
PUSH  REFRACTED  RAY  ON  STACK 
end  if 
end  else 
else 

if  (CURRENT  RAY  TYPE  =  VIEW  RAY)  then 
SET  INTENSITY  TO  BACKGROUND  COLOR 
end  else 

intii  STACK  EMPTY 
DISPLAY  PIXEL 

end  FOR-X 
end  FOR-Y 


Figure  5.2  -  Pseudocode  Description  of  the  Ray  Tracing  Process 


source  ray  is  on  the  bottom  of  the  ray  tuple  with  the  reflected  ray  above  it  and 
the  refracted  ray  above  that.  In  order  to  set  these  values,  it  is  necessary  to  pop 
the  stack  to  gain  access  to  the  source  ray.  Once  these  values  are  set  in  the  source 
ray,  the  program  returns  to  the  beginning  of  the  loop  and  pops  the  next  ray  off 
the  stack  discarding  the  last  ray  as  it  is  no  longer  needed.  [Ref.  1:  p.  377] 

The  intensity  algorithm  [Ref.  1:  p.  377]  although  very  simple  in  design, 
depending  on  the  number  of  light  sources  in  the  picture,  can  also  become  a  time 
intensive  part  of  the  ray  tracing  program.  The  entire  ray  data  set,  listed  in  Table 
3.1,  is  sent  into  this  routine  as  well  as  the  pointers  to  the  object  and  light  source 
list.  As  in  the  intersection  routine,  rays  are  generated  here.  In  Rogers,  they  are 
referred  to  as  shadow  feelers  and  the  same  term  is  used  in  this  study.  These 
shadow  feelers  are  those  rays  represented  as  the  vectors  from  the  point  of 
intersection  to  the  light  source,  see  Figure  5.3.  They  are  used  to  determine  the 
intensity  contributed  to  that  point  from  that  light  source.  Once  these  rays  are 
generated,  they  also  pass  through  the  intersection  routines  in  order  to  determine 
which  objects,  if  any,  the  light  rays  pass  through  en  route  to  the  intersection 
point.  The  first  test  that  must  be  done  is  to  determine  if  any  of  the  objects  passed 
through  are  opaque.  If  any  are  opaque,  then  no  iight  reaches  the  intersection 
point  from  "hat  light  source.  That  point  :s  then  considered  *o  he  lying  In  ueep 
shadow  with  respect  to  that  light  source.  If  none  of  the  objects  intersected  by  the 
light  ray  is  opaque,  then  the  light  intensity  needs  to  be  attenuated  accordingly  for 
each  of  them.  This  attenuation  entails  multiplying  the  intensity  at  each  point  by 


Figure  5.3 

Shadow  Feelers  [Ref.  3:  p.  64] 


the  transmission  coefficient  of  the  object.  This  process  then  needs  to  be  repeated 
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for  each  light  source.  A  running  total  of  the  intensities  is  maintained  to  be 
included  in  the  final  calculation.  The  calculation  takes  into  account  the  ambient 
light  and  the  light  that  comes  in  along  the  reflection  and  refraction  rays.  This 
process  then  produces  the  final  intensity.  If  the  input  ray  type  is  a  view  ray  then 
it  is  displayed.  If  it  is  a  reflected  ray,  it  becomes  the  Is  value  in  the  source  ray.  If 
it  is  a  refracted  ray,  it  becomes  the  I(  value  in  the  source  ray.  [Ref.  1:  pp.  376-377] 
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VI.  IMPLEMENTATION 


The  prototype  was  written  in  BORLAND’S  Turbo  Pascal  and  implemented 
on  an  IBM  AT  clone.  The  program  is  2500  lines  long  and  takes  two  hours  to 
generate  a  scene  of  200x200  pixels.  The  scenes  generated  on  the  AT  were  then 
displayed  using  the  RGB  monitor  of  a  Silicon  Graphics  IRIS  2400  graphics 
workstation. 

The  main  focus  of  this  study  was  to  develop  a  prototype  ray  tracer,  which  by 
itself  is  just  a  hidden  surface  removal  technique.  The  secondary  consideration  was 
to  integrate  an  illumination  model  into  the  ray  tracer.  Because  of  this  focus, 
more  time  was  spent  examining  the  ray  tracing  algorithm  than  any  of  the  global 
illumination  models  that  could  have  been  integrated  along  with  it. 

The  top  three  scenes  in  Figure  6.1  tested  the  ability  of  the  prototype  to 
perform  as  a  hidden  surface  remover.  The  program  proved  successful  in  this  area. 
For  these  scenes,  a  stub  was  used  in  place  of  the  illumination  model  and  each 
scene  was  lit  using  only  ambient  light.  From  left  to  right  the  scenes  show:  An 
unobstructed  view  of  the  ‘hree  objects,  described  in  a  later  section:  the  cube  and 
sphere  suspended  a'oove  the  door  out  with  the  cuDe  partially  blocking  the  sphere; 
the  cube  and  the  sphere  sunk  part  way  into  the  floor  with  the  cube  still  in  front.* 

‘These  were  the  only  scenes  generated  to  test  for  hidden  surface  removal.  The  remainder  of 
the  tests  were  done  trying  to  integrate  a  global  illumination  model. 


The  testing  of  the  global  illumination  model  integrated  into  the  ray  tracer 
unveiled  some  problems.  The  two  bottom  scenes  in  Figure  6.1  are  representative 
of  the  successful  results.  The  first  problem  discovered  was  that  attenuating  the 
intensity  by  the  full  distance  between  the  origin  of  the  view  ray  and  its 
intersection  point  produced  totally  blackened  objects.  The  results  shown  in 
Figure  6.1,  the  bottom  row,  were  obtained  by  either  dividing  the  intensity  by  two, 
see  the  scene  on  the  left,  or  by  not  dividing  it  at  all,  see  the  scene  on  the  bottom 
right  of  the  object.  The  second  problem  can  be  clearly  seen  by  the  black  line  that 
runs  up  the  center  of  the  floor  on  the  bottom  right  scene.  This  result  along  with 
those  test  that  generated  shadows,  not  shown,  indicated  that  the  intensities  for 
the  floor  were  reversed.  The  black  line,  clearly  seen  in  color  Figure  5  and  just 
vaguely  visible  in  color  Figure  4,  is  actually  specular  reflection  and  should  be 
much  brighter  than  the  rest  of  the  floor.  In  those  scenes  where  shadows  were 
generated  the  shadows  were  also  brighter  than  the  rest  of  the  floor— just  the 
reverse  of  what  it  should  have  been. 

A.  INPUT 

The  test  data  used  in  this  study  produced  the  scenes  shown  in  Figure  6.1. 
This  "esr,  data  was  in  the  form  of  a  sequential  die  with  ‘he  data  structure  ouriinea 
in  Figure  2.2.  The  data  had  in  it  one  picture  record,  one  light  in  the  lights  array, 
and  three  objects  in  the  objects  array.  The  first  object  is  a  cube  which  contains 
one  subobject.  This  subobject  has  one  common  part.  The  common  parts  record 


e  Generated  Scene 


then  points  to  the  six  polygon  records  that  where  used  to  construct  the  cube. 
Each  one  of  these  records  has  its  own  vertex  arrays.  The  second  object  is  a 
sphere.  This  object  has  one  subobject,  and  the  subobject  has  one  common  part. 
Since  a  sphere  can  be  its  own  bounding  volume  and  since  it  is  not  constructed  of 
polygons,  then  the  object  chain  for  the  sphere  need  go  no  further  than  the 
common  parts  array.  The  third  object  is  the  floor.*  This  consists  of  one 
subobject,  one  common  part  and  one  polygon.  All  objects  in  the  scene  are  opaque 
and  have  a  highly  reflective  surface. 

B.  OUTPUT 

The  output  generated  by  the  ray  tracer  is  in  the  form  of  a  bitmap,  with 
values  for  each  of  the  red,  green,  and  blue  components.  These  values  range  from 
0  to  1.  To  display  these  on  the  RGB  monitor  of  the  IRIS,  each  red,  green,  and 
blue  component  is  then  multiplied  by  255  and  assigned  an  index  in  the  color 


VII.  CONCLUSIONS 


A.  AREAS  OF  FUTURE  RESEARCH 

The  ray  tracer  is  a  powerful  tool  in  computer  graphics.  In  its  original  design, 
it  produced  the  finest  rendered  pictures  at  that  time.  Since  then  there  have  been 
numerous  extensions,  among  the  most  widely  known  are  the  ones  by  Phong,  Blinn 
and  Newell.  Kay.  and  Whitted  each  of  which  have  further  enhanced  the 
performance  of  the  ray  tracer  [Ref.  2:  pp.  343-344].  There  are  two  main  areas  for 
future  research:  global  illumination  models,  and  intersection  algorithms.  Both  of 
the  areas  are  of  great  interest  in  the  graphics  world.  Since  ray  tracers  and  global 
illumination  models  can  be  integrated  so  easily,  working  on  either  problem  would 
undoubtedly  lead  to  insight  into  the  other.  This  would  also  be  very  easy  to  do 
because  ray  tracing  naturally  lends  itself  to  a  modular  design  making  it  easy  to 
establish  hooks  for  the  testing  of  a  large  number  of  algorithms,  both  illumination 
and  intersection. 

B.  CONCLUSIONS 

■Vp  :iave  pxamineu  "fie  ‘nree  major  ireas  u  -av  'racing:  p-ene  :ara 

needed,  the  intersection  problem,  and  the  intensity  problem.  The  data  structure 
used  was  adapted  from  [Ref.  3]  and  proved  to  be  useful.  The  intersection 
problem,  although  involved,  is  not  complex.  The  algorithms  used  in  this 


implementation  are  simple.  The  inclusion  of  a  simple  global  illumination  model, 
was  easy  to  integrate  and  provided  fair  results.  The  ray  tracer  provides  an 
excellent  test  bed  program  and  implemented  can  provide  a  useful  tool  to  study 
numerous  problems,  not  only  in  lighting  and  shading  but  also  in  intersection 
determination. 
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APPENDIX  A  -  SOURCE  LISTINGS 


DECLARATIONS 


{  THIS  IS  THE  CONSTANTS  DECLARATION  SECTION  } 
maximum  size_of_stack  -  100; 
initial  pixel  z  =  0; 


{  TYPES  DECLARATION  AREA  } 

verticesarray  =  array  [  1. .4]  of  real; 

polygon  rec  =  record 
numvertices  :  integer; 

vertice  x,  vertice  y,  Venice  z  :  venices_array; 
surface_normal  x,  surface  normal  y,  surface_normal_z  :  real 
end; 


polygon_ptr  =  ‘polygon_array; 
polygon_array  =  array  jl..6]  of  polygon _rec; 


commonpartjrec  =  record 
K  ar,  K_ag,  K  ab  :  real;  {  AMBIENT  DIFFUSE  COEFFICIENT  } 

K  dr,  K_dg,  K_db  :  real;  {  DIRECT  DIFFUSE  COEFFICIENT  } 

K_sr,  K_sg,  K_sb  :  real;  {  SPECULAR  COEFFICIENT  } 

K  tr,  K  tg,  K  tb  :  real;  {  TRANSMISSION  COEFFICIENT  } 
obj  refraction  index  :  real;  {  OBJECTS  REFRACTION  COEFFICIENT  } 
obj  phong  exp  :  integer;  {  PHONG’S  SPECULAR  EXPONENT  } 
num_polygons  :  integer; 
polygons  :  polygon_ptr 
end; 

common_part_ptr  =  ‘common  part  array; 

common  part  array  =  array  (l..3|  of  common  part  rec; 


sub  object  rec  =  record 


rum  "ommon  Darts 


ommon  Darts 


:omrnun  Dart  ur: 


juo  bsphere  radius:  real;  ,  RADII  S  OK  ~1  BOH.JE<T‘  BOI  M.'l' 


•  ’HI  'I 


sub  bsphere  x 
sub_bsphere_y 
sub  bsphere  z 
subobj  type 
end; 


:  real.  {  CENTER  OF  BOUNDING  SPHERE  } 

:  real; 

:  real, 

:  integer;  {  0  SPHERE  1  PLANAR-POLYGON 


sub  object  ptr  —  “sub  object  array, 


1  'j 


»  ,Sj  O  »•  .-/V/V',- V 


sub  object  array  =  array  1  3  of  sub  object  rec 

object  rec  =  record 

num  sub  objects  integer 
sub  objects  sub  object  ptr. 

obi'bsohere  radius  real  <  RADII'S  OF  ORJECTS  ROUNDING  SPHERE  1 
ooj  nspnere  c  real.  ,  1  _'EN  7ER  >E  BU’1  N’>|\G  -i  ‘IKHh 

obj  bsphere  y  real, 

obj  bsphere  z  real, 

opcode  integer  {  CURRENTLY  NOT  USED  } 

end; 

object  ptr  =  ‘object  array 

object  array  =  array  1  4  of  object  rec. 

light  rec  =  record 

I  r  I  g,  I  b  real,  {  INTENSITY  OF  THE  LIGHT  ) 

light  x.  light  y,  light  z  real  {  LIGHT  POSITION  } 

dimension  1  dimension^  '■•a,  1  SOT  I  'hi'  ’■ 

end, 

light  ptr  =  ‘light  array. 

light  array  -  array  1  3  of  light  rec 

picture  rec  =  record 

num_objS  integer 

objects  object  ptr 

num  lights  integer 

lights  light  ptr 

global  refraction  index  real 

no  zero  real, 

ambient  r  real, 

ambient  g  real 

ambient  b  real 

background  color  r  real 
background  color  g  real 
background  color  b  real 
view  position  x  real 

view  position  y  real 

view  position  z  real 

screen  max  x  integer 

n>  e?*r 


'twi  max 


rav  rec 


rword 


ray  ivp«* 

rayty  pe 

ray  origin  x 

real 

{  >RIGIN  OF  RAY  } 

ray  origin  y 

real. 

ray  origin  i 

r»al. 

ray  v^rior  x 

'»al. 

{  DIRECTION  OF  R  AY  } 

i  >.  •  r 

'ii 

ray  vector  j 

real 

ray  slype 

ray  ly  pe 

|  TYPE  OF  SOI  RCE  RAY  } 

inieraertion  flag 

boolean 

obj  ldx 

inleger 

{  PATH  DECSR1B1NG  OBJECT  INTERSECTED  } 

a  u  S'  1  h  1  1  d  X 

;  n  t 't'r 

■•part  ldx 

integer 

polygon  idx 

integer 

.nteraertiun  x 

real 

{  INTERSECTION  POINT  ) 

intrra^''!  mn  y 

real 

inierxertion  t 

real 

DIST  AV’E  BETWEEN  RAV  S  ORIGIN  AND  INTERSECTION  POINT  } 


INTENSITY  <>F  light  coming  in  along  the  refracted  ray  generated  by  this 

HO 


i  ;r  i  Ig  i  to  r*ai 

INTENSITY  OI-  LIGHT  COMING  IN  ALONG  THE  REFLECTED  RAY  GENERATED  BY  THIS 
R  O 


I  sr  I  *g  I  sb  rrai 
r a>  .ink  ray  pi  r 

»nd  '  ray  pir  i 

•  P*i 

. \  AH . | 

v  ar 


njt  ti  ir  ini 
sy  sin  l<*»t 

IVK  III  !  *  X  I 

I  -HGNI  ON VERT ING  ALL  RAV  S  VECTOR?'  TO  I  NIT  VECTORS  } 

x  s  i  r*a! 


'ii  n  ii  m  r>ri  ni'-g'-r 

I  -ED  i  H<  -H  I  Hand  he<  AI  'E  THE  OHJEGT  P\TH>  GET  LONG  } 

pari  pai  H  or1  mon  pari  p'r 
»ijt«,r_  pal  t  nit  tij*-'  i  pi  r 

■  'il'  ir  '  •  < ’T i  y  p»- 
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{  USED  TO  IMPLEMENT  THE  STACK  } 
ray  _top, 
rav  _current, 
ray  next  :  ray  ptr; 


intensity  red. 


miensitv  >reen. 


intensity  _blue  :  real. 


pixel  _x, 
pixel  y, 

pixel  2  integer: 


ray  generation  number  :  integer; 


intersection  _point_x, 
intersection  point_y, 
intersection  point  i  :  real; 


old  intersection  point  x, 
old  intersection  point  _y, 
old  intersection  point  i  :  real; 


temp  integerl, 
temp  integer2. 
temp  integers  :  integer; 


reflected  ray  x, 
reflected  ray  y, 
reflected  ray  i  :  real, 


refracted  _ray_x, 
refracted  ray  _y. 


refracted  ray  *  : 


surface  normal  x. 
surface  normal  y, 
surface  normal  t  :  real; 


{  USED  TO  INITIALIZE  RAYS  } 

initial  ray  type  :  raytype; 

'mtial  rav  origin  x.  initial  rav  origin  v  initial  rav  origin  2  real: 


mum  -xv  *•••1.  or 


mum  *xv  •  ec'or 


nmai  xv  vector 


r,  it  i, xi  -xv  i  \  pe  -avr  v rje 

initial  intersection  tiag  oooiean. 


initial  obj  idx 
initial  subobjjdx 
initial  cpart  idx 
initial  polygon  idx 
initial  intersection _x, 
initial  intersection  y, 
initial  intersection  2 


integer; 

:  integer; 

:  integer; 

:  integer; 


* 

N 

»: 

i! 

N 


* 


\ 

% 


■« 

H 

H 

I) 

ll 


•l 

1 

«s 

*« 

'4 

i 

H 


initiald  :  real; 

initial  I  tr,  initiall  tg,  initial  I  tb  :  real; 
initial_I  sr,  initial_I_sg,  initial  l  sb  :  real; 

{  HOLDS  CURRENT  RAY  WHILE  IN  ACTUAL  RAY  TRACING  LOOP  } 
current  _ray  _type  :  raytype: 

current_ray_origtn_x,  current_ray _origin  y,  current  ray  _  jrigin_t  .  real; 

current  ray  vector  x,  current  ray  vector _y,  current  ray  vector  *  :  real; 

current  ray_stype  :  raytype; 

current  intersec tion_flag  :  boolean; 

current  objidx  :  integer; 

current  subobj  idx  :  integer; 

current  cpart  idx  :  integer; 

current  polygon  idx  :  integer; 

current  intersection  _x, 

current  in tersection  y, 

current  intersection  *  :  real; 

current  _d  :  real; 

current  I  tr.  current  I  tg,  current  I  tb  real; 
currentlsr,  current  l  sg,  current  l  sb  :  real; 

reflected  ray  :  boolean; 
refracted  ray  :  boolean; 

source  ray  num  ;  integer; 
sourceraytype  ;  raytype; 

{  POINTER  TO  PICTURE  RECORD  } 

picture  :  picture  rec; 

{ 

USED  TO  ESTABLISH  LINKED  LIST  OF  OBJ,  SUBOBJ,  LIGHTS,  CPARTS,  AND 
POLYGONS. 

} 

light  cntr  ;  integer; 
lightcurrent  ;  lightptr; 

obj  cntr  :  integer; 
objcurr  :  object  ptr; 


<• 

A 

A 


subobjcntr  :  integer; 
subobj  curr  •  sub  object  Dtr: 

part  rur  integer; 
cpart  curr  :  common  part  ptr; 


polycntr  ;  integer; 
poly  next, 

poly  curr  :  polygon  ptr; 


vertice  cntr  :  integer; 


{  USED  WHEN  READING  IN  DATA  } 
poly  loop _cnt  .  integer; 
light  loopcnt  :  integer; 
cpnrtloopcnt  :  integer; 
subobjloopcnt  :  integer; 
object  loop  cnt  :  integer; 
vertice  ioop  cnt  integer: 

intersection _X,  intersection  Y,  intersection _Z  :  resl; 

{  USED  TO  IDENTIFY  INTERSECTED  OBJECTS  } 
object  idx. 
subobj  idx, 
cpnrt  jdx, 

polygon  idx  integer; 


intersection  JUg  ;  boolesn, 


INTERSECTION  PROCEDURES 


i 


{ 


. * . INTPRCSC  PAS  * . * . * . 

*  These  a re  the  intersection  procedure*  used  in  the  ray  tracer  program 


\ 


. SPHERE  INTERSECTION . . . 

*  CALLED  FORM  CHECK  FOR  SUBOBJ- INTERSECTION 

*  r  XLLS  TO  NONE 

*  DESC  CALCULATES  THE  INTERSECTION  POINT  BETWEEN  A  RAY  /VECTOR  AND  A 

*  SPHERE 

*  INPUT  The  centerpoint  of  the  sphere  and  it's  radius  The  direction  of 

*  the  ray  and  a  known  point  on  the  ray  --  which  would  be  it's 

*  origin 

*  OUTPUT  A  flag  indicating  whether  or  not  there  was  an  intersection  and 


ner*  was  'he  in  usi  nterser'ion  ;i«>int  u«elf 


I 


c  :=  (  sqr(Cx)  +  sqr(Cy)  +  sqr(Ci)  + 
sqr(Px)  +  sqr(I^)  +  sqr(Px)  - 

(2  *  Cx  *  Px)  -  (2  *  Cy  *  Py)  -  (2  *  Cx  *  Px)  -  sqr(r)); 
radical  :=  (tqrt  (  sqr(b)  -  (4  *  a  *  c))); 

(  START  COMPUTATIONS  ON  QUADRATIC  EQUATION  } 

if  (radical  <  0)  then  begin 

writeln(’ WARNING  -  un  agio  nary  number  possible  in  SPHERE  INTERSECTION’); 
o  intersection  flag  :=  false 
end 

else  begin 

if  (radical  =  0)  then  begin 

{ 

IF  0  THEN  JUST  ONE  INTERSECTION  POINT(THE  LINE  IS  TANGENT  TO  THE  SPHERE) 

} 

tl  :=  (-b  /  2  •  a);  {  SOLVE  FOR  t  } 

Sphere  x  :=  Px  +  (Vx  *  tl);  {  CALCULATE  POINT  USING  t  } 

Sphere  y  =  Pv  -  (Vy  *  tl). 

Sphere  s  :=  Px  -+•  iVx  *  tl), 

©intersection  flag  :  =  true 
end 

else  begin 

{ 

THERE  WERE  TWO  INTERSECTION  POINTS  -  ONE  ENTERANCE  POINT  AND  ONE  EXIT 
POINT  SOLVE  FOR  BOTH  ts’ 

} 

tl  -  (-b  -  (sqrt  (  sqr(b)  -  (4  *  a  *  c))))  /  2  *  a; 
t2  .=  (-b  -  (sqrt  (  sqr(b)  -  (4  *  a  *  e))))  /  2  *  a; 

XI  :=  Px  -  (Vx  *  tl). 

Y 1  .=  Py  *  (Vy  *  tl); 

Zl  :=  P*  ♦  (Vx  *  tl), 


< 

CALCULATE  DISTANCE  FOR  BOTH  INTERSECTION  POINTS  FROM  THE  POINT  OF  ORIGIN 

I 

distancel  =  (tqrt  (sqr(Xl  -  Px)  +  sqr(Yl  -  Py)  +-  »qr(Zl  -  Px))); 

X2  --=  Px  -  (Vx  *  t2) 

Y2  =  Py  -  (Vy  *  t2) 

'2  -  Pi  -  'Vx  *  .21 

UMan.-eJ  j  rt  -or.  \2  i’xi  -qn  V2  i’"i  -onZ'J  ’ill! 

{  COMPARE  DISTANCES  AND  SELECT  THE  INTERSECTION  POINT  THAT  IS  CLOSER  ) 
if  distancel  -  distance2  then  begin 
Sphere  x  -  X 1 . 

Sphere  y  =  Y  1 , 

Sphere  s  -  Zl, 
end 

else  begin 
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Sphere_x  :=  X2; 
Sphere  _y  :=  Y2; 


Sphere  _*  :=  Z2; 
end; 

o  intersection  flag  :=  true 

end:  <  ELSE  } 

diffx  :=  Px  -  Sphere  x; 
diffy  :=  Py  -  Sphere  y; 
diff*  :=  Pj  -  Sphere  *; 

{ 

A  CHECK  TO  INSURE  THAT  THE  POINT  SELECTED  ISN’T  THE  RAYS  POINT  OF  ORIGIN 

} 

if  ((diffx  <=  0.0000)  nnd  (diffy  <=  0.0000)  and  (diff*  <=  0.0000))  then 
o  intersection  flag  :=  false; 

end  {  ELSE  ) 

end;  t*  SPHERE  INTERSECTIONS  *} 


(PA) 

1 

. .  CALCULATE  PLANE  EQUATION  * . 

*  CALLED  FROM  FIND  INTERSECTED  POLYGON 

*  CALLS  TO  NONE 

*  DESC  Calculates  the  constants  of  the  equation  of  a  plane  when  given 

*  three  point  on  the  plane 

*  INPUT  Three  vertices  of  a  planar  polygon 

*  OUTPUT  The  A,  B,  C,  and  D  constants  for  the  equation  of  a  plane 


) 

procedure  calculate  plane  equation(Xl,  Yl,  Zl  real; 

X2,  Y2,  Z2  real 
X3,  Y3,  Z3  real, 
var  A,B,C,D  real  ), 


begin 


4  Yl  * 

IZ2 

731  •  Y2 

*  (Z3  -  Zl) 

-  Y3  *  (Zl  -  Z21. 

i  ' ; 

\.’ 

\Z)  - 

■  \.:  ■  -iii 

•  "  Yi 

\i  • 

\  J 

N.r  ■  v: 

• 

i  \ ;;  •  \ 

'  H> 

t  •  9  •*  \ 

-  •  \  l 

•»  **•) )  ■  i 

\  J  \  \  /.if- 

■  (XS  *  II 

Yl  * 

72)  -  (Y2  * 

Zl))). 

»nd.  (  calculate  plane  equation  } 


**•«•*••***•****•*  FIND  INTERSECTION  point  ********»***»»«***»**«* 

CALLED  FROM:  FIND  INTERSECTED  POLYGON, 
CHECKFORSUBOBJINTERSECTION 
CHECKFORINTERSECTION 

CALLS  TO  :  NONE 

DESC  :  Calculates  the  intersection  point  between  a  ray /vector  and 
a  plane. 

INPUT  :  Ray  direction  and  a  known  point  on  the  ray,  i.e.,  it’s 
point  of  origin.  The  constants  (A,B,C,D)  of  the  equation 
of  a  plane. 

OUTPUT  :  The  intersection  point. 


procedur-  nnd_intersection_point(i_A,  i_B,  i_C,  i_D  :  real; 

irayx,  i_ray  y,  i_ray_s  :  real; 
i  source  x,  i  source_y,  i_source_s  :  real; 
var  o intersection  point_x, 
o  intersection  point  _y , 
o  intersection  point  z  real): 


t  .  real; 


{  SET  UP  FOR  FINDING  t  FROM  THE  EQUATION  FOR  A  LINE  AND  A  PLANE  } 
t  :=  (i_D  -  ((i_A  *  i  source  x)  +  (i_B  *  i  source  y)  +  (i_C  *  isources)))  / 

( (i _ A  *  i  ray  x)  + 

(i_B  *  irayy)  + 

(i  C  *  i_ray_s)); 

{ 

SUBSTITUTE  t  BACK  INTO  THE  EQUATION  FOR  A  LINE  TO  GET  THE  INTERSECTION 
POINT 

} 

o  intersection  point _x  :=  (i  ray  *  *  t)  +  i  source  x; 

o  intersection  point  y  :=  (i  ray  y  *  t)  +  i  source _y; 

o  intersection  point  i  :=  (i  ray  »  *  t)  +  i  source  s; 


end;  {  find  intersection  point  } 


CALLED  FROM:  FIND  INTERSECTED  POLYGON 
CHECK  FOR  SUBOBJ  INTERSECTION 
CHECK  FOR  INTERSECTION 
CALLS  TO  :  NONE 

DESC  ;  This  calculates  the  constants  A.B.C.D  of  the  equation  of  a  plane 
given  only  a  point  on  the  plane  and  the  surface  normal  of  the 
plane 

INPUT  A  point  and  the  surface  normal  of  the  plane  whose  equation  you 


*  are  trying  to  figure  out. 

*  OUTPUT:  The  constants  A,B,C,D  of  the  equation  of  a  plane. 


} 

procedure  ealeq(bsphere_x,  bsphere_y,  bsphere  i, 
vector  x,  vector  y,  vector  _*  :  real; 

var  B.  0.  D  :  real); 

begin 

A  :=  vectorx; 

B  =  vector  y: 

C  :=  vector  i; 

D  :=  (vector_x  *  bsphere_x)  - 
(vectory  *  bsphere  y)  •+• 

(vector  1  *  bsphere _*); 

end:  {  caleq  } 


{ 


. *****  POLYGON  ORIENTATION  **** . 

CALLED  FORM:  FIND  INTERSECTED  POLYGON- 
CALLS  TO  NONE 

DESC  This  checks  to  see  if  the  incoming  ray  will  hit  the  front  face  of 
this  polygon  It  does  this  by  comparing  the  angle  between  the 
surface  normal  of  the  polygon  and  the  inverse  of  the  incoming  ray 
If  the  angle  is  greater  than  90  degrees  then  the  ray  is 
approaching  the  back  of  the  polygon 
INPUT  The  direction  of  the  incoming  ray  The  surface  normal  of  hte 
object 

OUTPUT  A  boolean  value  TRUE  FALSE  depending  on  whether  or  not  the 
polygon  is  facing  the  right  direction 


procedure  polygon  orientation  (view  vector  x  view  vector  >  view  vector  <  real 
surface  normal  x,  surface  normal  y.  surface  normal  i  real, 
var  o  good  orientation  boolean) 

var 

<-i-wine  theta  real 

ft  .  >r  **  *| 

-•  -  f  ft  •  i r*  %  •*  »■  r*n  *  --\i 

:<  ••  ^r**  iu*  ’ 

begin 

< 

r  \hi  nil  dot  pkodi  <  t  oi  rm  i\\  frsi  oi  i  hf  ml  vn  \  m  mu  \m»  thi  d  him  i 

NORMAL  OF  THE  PO|  YOON  IN  gl  ESTION 
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dot  product  :=  (-  view  vector_x)  *  surface  _normal  x 
(-  view  vector  _y)  *  surface  normal  y  -+- 
(-  view_vector  i)  *  surface  normal  z; 

{  CALCULATE  THE  MAGNITUDE  OF  THE  VECTORS  } 

length  view  vector  :=  sqrt(sqr(-  view  vector  x)  + 
sqr\-  view  vector  y)  — 
sqr(*  view_vector_z)); 

lengthsurfacenormal  :=  sqrt(sqr(  surface  normal  x)  ■+• 
sqr(  surface  normal  y)  + 

sqrf  surface_normal_z)); 

{  CALCULATE  THE  COSINE  OF  THE  ANGLE  BETWEEN  THE  RAYS  } 

cosine_theta  :=  dot_product  /  (length_view_vector  *  length _surface_normaJ) 

if  (cosine  theta  >  0)  then 
o  good  orientation  :  -  true 
else 

o  good  orientation  :=  false; 
end.  {  polygon  orientation  } 


{  PA} 

{ 

. FIND  INTERSECTED  POLYGON . 

*  CALLED  FROM  CHECK  FOR  SCBOBJ  INTERSECTION 

*  CALLS  TO  POLYGON  ORIENTATION 

*  CALEQ 

*  FIND  INTERSECTION  POINT 

*  CALCULATE  PLANE  EQUATION 

*  DESC  This  determines  if  there  is  an  intersection  between  •  line  ray 

*  and  a  polygon  If  there  is  it  calculates  what  it  is 

*  INPUT  The  direction  of  the  shooting  ray.  a  known  point  on  that  -a> 

*  it  s  origin  and  the  object  path  identifying  the  subobiert  to 

*  examine 

*  Ol  TP'  T  A  flag  indicating  whether  or  not  a  polygon  was  hit  If  >ne  w  %< 

hit  the  path  identifying  which  >ue  it  was  and  the  actual 

*  intersection  point  itself 


.  •  .  iv  ,  , 

i  •  -  i  r 

t>  'll 

niV’hi  ids  integer 
'  er  -  f~lvg  -n  -ntersec  t  i..r  « 
polsg.m  .nlersection  \ 
poitg  >n  ,nlersec|  .on  t  -»« 

'  ar  par*  id  s 

[» .l>  «  >n  ids  ini  eg  r* 

>  ar  n terse.- 1  t..p  flag 


*V 


cpart  cnt  •  =  1; 
polygon  cnt  :=  1; 
vertice  cnt  :=  1; 
intersection  _found  :=  false; 
good  orientation  :=  false; 

cpart  path  :=  picture  objects  *|i  obj  idxl.sub  objects'll  subobj  idx|. 
oramon  parts: 

polygon  path  :=  cpart  path '[cpart_cnti. polygons; 

( 

THIS  CHECKS  EACH  COMMON  PART  FOR  AN  INTERSECTION  WITH  ONE  OF  IT’S 
POLYGONS 

} 

repeat 

\  THIS  LOOP  CHECKS  EACH  POLYGON  OF  A  COMMON  PART  FOR  INTERSECTION  } 

repeat 

point  outside  polygon  :  =  false; 

Mvivgon  orientationli  ray  x.  i  ray  y.  i  ray  i, 

polygonpath  jpolygon_cntj.surtace_normai_x, 
polygon _path  ‘[polygon  cnt j. surf ace_normal_y, 
polygon  path  '[polygon_cnt].surface_normal  z, 
good  orientation); 

J  good  orientation  then  begin 

caleq( polygon  path  '  polygon  cnt j.vertice_x[l], 
polygon  path  '(polygon  cnt]. vertice_y|l], 
polygon  path  '[polygon_cnt].vertice  z[l], 
polygon  path  '[polygon_cnt].surface_nornial_x, 
polygon  path  ' [polygon  cnt]. surface_normal_y, 
polygon  path  '!polygon_cnt].surface_normal  z, 

A.B.C.D); 

find  intersection  point(A,  B,  C,  D, 
i_ray_x, 
irayy, 
i_ray_z, 

i_source_x,  i_source  y,  i  source  z, 
polygonX,  polygonY,  polygonZ); 

HKCK  r'r  MAKE  SI  RE  YOi;  \RE  NOT  CONSIDERING  THE  SOURCE  POINT(ORIGIN)  OF 
i  HE  R.\  Y 

) 

diffx  :=  i_source_x  -  polygonX; 
diffy  :=  i  source  y  -  polygonY; 
diffz  :=  i  source  z  -  polygonZ; 

{ 

THIS  SETS  A  FLAG  IF  YOU  DO  CONSIDER  THE  SAME  POINT,  THE  CHECKS  IN  THE 
OTHER  PROCEDURE  SHOULD  PREVENT  THIS  BUT  JUST  IN  CASE. 


} 


if  ((diffx  <=  0.000)  and  (diffy  <=  0.000)  and  (diffi  <  =  0  000)) 

then  begin 


't*VViVr 


{  writeln(sysout,’set  trip-wire’);} 
intersections!  l!  :=  false 
end 

else  begin 


THIS  LOOP  CHECKS  EACH  EDGE  OF  THE  POLYGON  TO  SEE  IF  THE  INTERSECTION 
POINT  LIES  INSIDE  OR  OUTSIDE  OF  IT. 


repeat 


{  THIS  SELECTS  THE  FIRST  VERTEX  OF  A  POLYGON  } 

xl  :=  polygon_path  *  [polygon  cnt].vertice_x[vertice_cntj; 
yl  :=  polygon_path  “(polygon_cntj.vertice_y[vertice_cnt]; 
si  :=  polygon_path  “[polygon_cnt], vertice  z[vertice  cntj; 


if  (vertice  cnt  =  (polygon  paoh  '  polygon  cnt;.num  vertices  - 

1)) 


WHEN  YOU  PICK  THE  NEXT  TO  LAST  VERTEX  YOU  CAN  SELECT  THE  NEXT 
CONSECUTIVE  VERTEX  TO  ESTABLISH  THE  EDGE  THROUGH  WHICH  YOU  WANT  THE 
BOUNDING  PLANE  TO  PASS.  YOU  THEN  MUST  PICK  ONE  OF  THE  OTHER  VERTICES  (AND 
IT  DOES  NOT  MAKE  ANY  DIFFERENCE  WHICH  ONE.  I  PICK  THE  FIRST  ONE.)  TO  BE 
USED  TO  PUT  INTO  THE  EQUATION  OF  THE  PLANE  THE  RESULT  OF  WHICH  IS  COMPARED 
AGAINST  THE  RESULT  THAT  COMES  FROM  PLUGGING  THE  INTERSECTION  POINT  INTO 
THE  EQUATION  OF  THE  BOUNDING  PLANE. 


then  begin 

x2  :=  polygon  path  "[ polygon  cntj. 

vertice_x[vertice_cnt  +  1 J; 
y2  :=  polygon_path  “[polygon  cnt], 
vertice_y[vertice_cnt  +  l); 
z2  :=  polygon_path  “(polygon_cnt], 
vertice  z [vertice  cnt  +  1  ] ; 


x3  :=  polygon  path  *[polygon_cnt].vertice_x[ij; 
y3  :=  polygon_path ‘[polygon_cnt|. vertice  y[lj; 
z3  :=  polygon_path '|polygon_cntj.vertice_z[lj 
end 

else  begin 


IF  YOU  DO  NOT  HAVE  THE  NEXT  TO  LAST  EDGE  THEN  .JUST  SELECT  THE  NEXT 
CONSECUTIVE  VERTICE  TO  ESTABLISH  THE  EDGE  FOR  THE  BOUNDING  PLANE  A  ND 
THE  ONE  AFTER  THAT  TO  PLUG  INTO  THE  EQUATON  OF  THE  PLANE. 


x2  :=  polygon_path  ‘(polygon_cntj. 

vertice_x[vertice_cnt  +  l] ; 
y2  :=  polygon _path  “(polygon _cnt|. 

vertice_y[vertice_cnt  +  1|; 
z2  :=  polygon  path  *jpolygon_cnt|. 


vertice  i  vertice  cnt  *  l 

xS  :=  polygon  path  '  polygon  cnt 
vertice  _xl  vertice  cm  »  2 
y3  -  polygon  path  '  polygon  cnt 
vertice  yivertice  cnt  ♦  2 
1$  -  polygon  patn  poiyg'  n  nt 

vert  ice  *  vert  ice  cnt  -  2 

end; 

if  (vertice_cnt  =  polygon _path  '  polygon  cnt  num  verticet  ) 

{ 

IF  THE  VERTICE  SELECTED  IS  THE  LAST  ONE  THEN  JUST  PICK  THE  FIRST  VERTJCE 
TO  ESTABLISH  YOUR  BOUNDING  EDGE  AND  THE  SECOND  VERTICE  TO  PLUG  INTO  THE 
EQUATION  OF  THE  PLANE 
} 

then  begin 

x2  :=  polygon_path  '  polygon  cnt  vertice  x  1  , 
y2  =  polygon  path"  polygon  ;ni  vertice  y  1 

z2  :=  polygon  path  polygon  cnt  vertice  i  1 

x3  :=  polygon  path  ‘  polygon  cnt  vertice  x  2  . 

y3  :=  polygon  path '  polygon  cnt  vertice  y  2  . 

z3  :=  poiygon_path  polygon  cnt  vertice  t  2 
end; 


{ 

ESTABLISH  THE  ARBITRARY  POINT  THROUGH  WHICH  THE  PLANE  WILL  PASS 

} 

anchor  x  :=  ((x2  -►  xl)  /  2)  10; 
anchor  y  :=  ((y2  -+•  yl)  /  2)  +  10; 
anchor  *  :=  ((*2  +  il)  /  2)  10; 

calculate  plane _equation(  xl,  yl,  zl, 
x2,  y2,  i2, 

anchor_x,  anchor_y,  anchor  i, 

A,  B,  C,  D); 

{ 

THE  RESULT  OF  PLUGGLNG  IN  THE  VERTICE  OF  THE  POLYGON  INTO  THE  EQUATION 
OF  THE  PLANE. 

) 

•narkerD  -  A  *  x«l  *  13  *  ..»!  -  C  *  ;31 

THE  RESULT  OF  PLUGGING  IN  THE  INTERSECTION  POINT  INTO  THE  EQUATION  UK 
THE  PLANE. 

} 

checkpointD  :=  (A  *  polygonX)  + 

(B  *  polygonY)  + 

(C  *  polygonZ); 

if  ((markerD  <=  -D)  and  (check  point  D  <  =  - D) )  then 
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IF  THE  RESULTS  HAVE  THE  SAME  SIGN  THEN  THEY  BOTH  LIE  ON  THE  SAME  SIDE  OF 
THE  BOUNDING  PLANE  HENCE  THE  INTERSECTION  POINT  LIES  WITHIN  THE  POLYGON 
W  ITH  RESPECT  TO  THAT  EDGE 


intonations  vertice  cnt 


true 


I 

IE  THE  RESULTS  DON'T  HAVE  THE  SAME  SIGN  THEN  THEY  LIE  ON  OPPOSITE  SIDES 
OF  THE  BOUNDING  PLANE  AT  THIS  POINT  THE  INTERSECTION  POINT  HAS  BEEN 
PROVEN  TO  LIE  OUTSIDE  THE  POLYGON 

if  ((markerD  >  =  -D)  and  (checkpoint  D  >—  - D ) )  then 
intersection*  verticecnt,  .  =  true 
else 

intersections! vertice  cnC  :=  false; 
vertice  cnt  —  vertice  cnt  1 

in'u  i  vert  in*  m  polygon  path'  polygon  cnt  num  vertices), 
end 
end. 


CHECK  THE  POLYGON  INTERSECTION  ARRAY  TO  SEE  IF  THE  INTERSECTION  POINT 
FAILED  THE  INSIDE  TEST  FOR  ANY  OF  THE  EDGES. 

> 

for  loop  cnt  ■=  1  to  polygon  path ‘'polygon  cnt,  num  vertices  do 
if  not(intersections[loop  cnt|)  then 
point  outside  polygon  =  true; 

if  point  outsidepolygon  then 
intersection  found  :=  false 
else 

intersection  found  :=  true; 

polygoncnt  :=  polygon_cnt  +  I; 
until  ((polygon  cnt  >  cpart_path  *  cpart_cntj.num_polygons)  or 
(intersection  found)), 

cpart  cnt  :=  cpart  _cnt  f  1; 

until  ((cpart  cnt  >  picture. objects  ‘]i  obj_idxj.sub  objects  *|i  subobj  idx). 
lum  "ommon  parts)  or  intersection  found)): 

'ET  up  THE  OUTPUT  FOR  THE  PROCEDURE  1 

o  cpart  idx  :=  (cpart  cnt  -  1); 
o  polygon_idx  :=  (polygon  cnt  -  I); 


end,  {  find  intersected  _polygon  } 


{.PA} 

{ 

******* . *  CHECK  FOR  SUBOBJ  INTERSECTION  . . *** 

*  CALLED  FROM:  CHECK  FOR  INTERS ETION 

*  CALLS  TO  :  CALEQ 

*  FIND  INTERSECTION  POINT 

*  SPHEREINTERSECTIONPOINT 

*  FIND  INTERSECTIONPOINT 

*  DESC  :  Check  to  find  out  if  the  shooting  ray  intersects  this 

*  subobject’s  bounding  volume. 

*  INPUT  :  The  object  whose  bounding  volume  has  been  hit. 

*  The  direction  of  the  shooting  ray. 

*  The  origin  of  the  shooting  ray. 

*  OUTPUT  :  A  flag  indicating  whether  or  not  there  has  been  an  intersection. 

*  The  intersection  point  --  if  there  is  one. 

*  The  path  to  the  intersected  object. 

*••*•**••*»********•**•»*********•*••*»*******#******************»*»*< 

} 

procedure  check_for_subobj  intersection  (i_object_idx  :  integer; 

i  ray  x,  i_ray_y,  i  ray  z  :  real; 
isourcex, 
i  source  y, 
i  source  *  :  real; 
var  ©  intersection _x, 
ointersectionjy, 
o  intersection  i  ;  real; 
var  o  subobj  idx, 
o  cpart  idx, 
o  polygon  idx  :  integer; 
var  o_intersection_flag  :  boolean); 

var 

closest  _object  :  real; 

{  VECTOR  BETWEEN  THE  RAYS  ORIGIN  AND  INTERSECTION  POINT  } 
view  polygon_vector  x, 
view_polygon_vector_y, 
view  oolvgon  vector  i  real: 

iistance  from  intersection. 

distance_from_view_position  :  real; 

subobj  path  :  sub  object  ptr; 
obj_path  :  object  _ptr; 

cpart  cnt, 
polygoncnt, 

subobj  cnt  ;  integer;  {  USED  TO  GO  INTO  RESPECTIVE  ARRAYS  } 
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A,  B,  C,  D  :  real; 


{  CONSTANTS  FOR  EQUATION  OF  A  PLANE  } 


{ 

INTERSECTION  POINT  BETWEEN  THE  RAY  AND  PLANE  THE  BOUNDING  CIRCLE  IS 
INSCRIBED  ON 
} 

bplane  intersection 
bplane  intersection  y, 
bplane_intersection  z  ;  real; 

{ 

INTERSECTION  POINT  BETWEEN  THE  RAY  AND  A  POLYGON 

} 

polygonX,  polygonY,  polygonZ  :  real; 
intersection  flag  :  boolean; 

begin 

o  3ubobj_idx  :=  0; 
ocpartidx  ;=  0; 
opolygonidx  :=  0; 
ointersectionflag  :=  false; 
o  intersection _x  :=  0.0; 
ointersectiony  :=  0.0; 
ointersectionz  :=  0.0; 

{ 

ESTABLISH  A  DEFAULT  DISTANCE  WITH  WHICH  THE  ACTUAL  DISTANCES  WIL  BE 
COMPARED. 

} 

closest  object  :=  10000.0; 

subobj  path  :=  picture. objects *(i  object_idx].sub_objects; 
obj_path  :=  picture. objects; 

{  THIS  LOOP  CHECKS  EACH  SUBOBJECT  OF  AN  OBJECT  } 
for  subobj  cnt  :=  1  to  obj  path  *(i_object_idxj.num_sub_objects  do  begin 

{  FIRST  ESTABLISH  THE  PLANE  ON  WHICH  TO  DRAW  THE  BOUNDING  CIRCLE  } 

raleqfsubobj  path*'subobj  cnt!. sub  bsphere 
subob|  path  "subobj  cnti.sub  bsphere  y. 
subobj  path  "subobj  '-nc  .sub  bsphere  z. 
i  ray  x,  i_ray_y,  i_ra.y _z, 

A,  B,  C,  D); 

{  FIND  THE  INTERSECTION  POINT  ON  THAT  PLANE  } 
find  intersection  point(A,B>C?,D, 

i_ray_x,  i_ray_y,  i_ray_z, 
i  source  x,  i  source_y,  i_source_z, 
bplaneintersectionx, 


bplane  inl»rw<  ti.ni  v 
bplanr  intersection  Cl 


{ 

DETERMINE  THE  DISTANCE  BETWEEN  THE  INTERSECTION  POINT  AND  THE  (  ENTER 
POINT  OK  THE  BOUNDING  SPHEH h  CIRCLE 

i 

distance  from  intersection 

sqrt(tqr(  tubobj  path  '  tubobj  cnt  sub  btphere  «  - 
bplane  intersection  x|  *- 
•qrfsubobj  path  '  tubobj  cnt  tub  btphere  y  ■ 
bplane  intersection  !  • 
sqr(tubobj  path  '  tubobj  cnt  tub  btphere  i  - 
bplatte  intersection  i)  ). 

{  IF  THE  DISTANCE  IS  LESS  THAN  OR  EQUAL  TO  THEN  YOU  HAVE  AN  INTERSECTION  } 
if  (distance  from  intersection  <  - 

tubobj  path  ~  tubobj  cnt  tub  btphere  radius)  then  begin 
if  (tubobj  path  '  tubobj  cnt  tubobj  type  =  0)  then  begin 

{ 

IF  SUBOBJECT  TYPE  IS  A  SPHERE  THEN  YOU  NEED  TO  USE  THIS  INTERSECTION 
PROCEDURE 

} 

sphere  intersect  ion  (i  source  x,  i  source  y,  i  source  i, 
i  ray  x,  i  ray  y,  i  ray  i, 
tubobj  path ‘  tubobj  cnt  tub  btphere  x, 
tubobj  path  ‘  tubobj  cnt'  tub  btphere  y, 
tubobj  path ‘  tubobj  cnt  sub  btphere  t, 
subobj  path  ‘;subobj  cnt  .sub  bephere  radius, 
intersection  flag, 
polygonX,  polygonY,  polygonZ); 
cpart  cnt  :=  1; 
polygon  cnt  0; 

end 

else  begin 

{  IF  IT  ISN’T  A  SPHERE  THEN  USE  THESE  } 

find  intersected_polygon(i_ray_x,  i  ray  y,  i  ray  i, 
i_source_x,  i_source  y,  i  source  t, 
i_object_idx, 
subobj  cnt, 

polygonX.  polygon  Y  polygonZ. 

Dart  -nr. 
polygoncnt, 
intersection  flag); 

end;  {*  ELSE*} 
if  intersection  flag  then  begin 

{ 

IF  THERE  HAS  BEEN  AN  INTERSECTION  THE  ESTABLISH  THE  VECTOR  BETWEEN  THE 
ORIGIN  OF  THE  RAY  AND  THE  INTERSECTION  POINT, 


V  *  '  i  '  H  I  *  -■  ■  ■  r  % 

v  i»m>  (h  J  v  ff  •  ‘f»  >•»■•*  .  r  > 


%  -  •  i*  r »  *■ 


.  i-  u «  'i  \ 

Si  MJf  .  f  V  h  A 

n»»  p*.  X  g- >11  x *■  ■  r  i  i  vur>  «■  i 

\ .  •  .  '  '  •»»  'h-  i  1  :  >  i  1  '  U  ,  ■ 

lisian  *•  fr  m  1 1»»  x.(rt '  *nr‘  *  lew  pi.lygm  ‘e<  for  xl  • 

«|rlvif»  •  1  y  ■  •  •  n  ‘fit'll  >1  • 

•tiriview  polygon  vector  ill 

if  !i»t«m*  'V  .m  i  if»  position  . loses!  oltijei  I  f  hr ti  begin 

I 

<  <>\U'\fU  II  AGAIN"!  rHKtirHKH  llKU\a>  AND  sKI  KT  I  UK  CLOSEST  ONE 

I 

loses!  >i>je<  t  distance  from  \  lew  pi. ml  imi 
•  iuImiI.j  ill  x  xulmfij  <  nl 
■  .nr'  i«  ■  nr1  ' 

polygon  nix  polygon  mi 

■  m lersec  i urn  iu«  mi»»rs*«  r  i« »n  fUtt 
.1  inlfrwrtiun  x  poly  gun  X 

.  inlfrxfrtu  n  v  poixgonA 

intfrxfiiii.il  <  polygon/ 

f  rid 
f  ml 
f  ml 

f nil ,  {*  FOR  *} 

fn«l.  1  check  f<ir  nbjfct  intersection  ) 

I  I’M 
f 

. CHECK  FOR  INTERSECTION . 

*  CALLED  FROM  MAIN  AND  CALC!  LATE  INTENSITY 

*  CALLS  TO  CALEQ 

*  FIND  INTERSECTION  POINT 

*  CHECK  FOR  SUBOBJECT  INTERSECTION 

*  DESC  DETERMINES  IF  THERE  IS  AN  INTERSECTION  BETWEEN  THE  SHOOTING  RAY 

*  AND  THE  BOUNDING  VOLUME  OF  AN  OBJECT 

*  INPOT  Direction  of  the  xhootmg  rav 

)rigin  it  he  -.nooiing  rav 
inter  nt. o  .  h»*  inject  trrav 

OUTPUT  :  Flag  indicating  wnetner  or  not  there  was  an  intersection. 

*  If  there  was  an  intersection  then  then  the  intersection  point. 

*  The  path  to  the  intersected  polygon. 

} 

procedure  check  _for_intersection  (i  ray  x,  i  ray  y,  i  ray  z  :  real; 

i_source_x,  i_source_y,  i_source_z  :  real; 

88 


y-. 


t  picture  objects  object  ptr 
vv  o  intersection  x. 
o  intersection  y, 
o  intersection  s  reel, 
var  o  object  jdx, 
o  subobj  idx. 

)  'port  nix, 
o  poly  ldx  integer; 
var  o  intersection  flag  boolean), 


’  ar 

'loses!  'bje'-t  r«*.»l 

object  cnl, 
object  idx, 
subobj  cnt, 
cperi  cnt, 

polygon  cnl  integer, 

distance  from  viewposition, 
distance  from  intersection  real, 


INTERSECTION  POINT  BETWEEN  SHOOTING  RAY  AND  THE  PLANE  THE  BOUNDING  CIRCLE 
IS  INSCRIBED  ON 


bplane  intersection  x, 
bplane  intersection  y, 
bplane  intersection  i  real; 

view  bplane  vector  x, 
view  bplane  vector  y. 
view  bplane  vector _i  :  real; 


INTERSECTION  POINT  BETWEEN  THE  SHOOTING  RAY  AND  A  POLYGON. 
polygonX,  polygon V',  polygonZ  :  real; 


CONSTANTS  FOR  THE  EQUATION  OF  A  PLANE. 


v.  B. 


\ 

found 


I 


D  real; 

intersection  :  boolean; 


distance  ;  real; 


o  poly  ldx  :  =  0; 
o_inlerseclion  _x  :  =  0.0; 
o_intersection_y  :=■  0.0, 
o  intersection  _*  =  0.0; 
o  intersection  flag  =  false: 
loseni  inject  -  1U00U  0, 

found  intersection  :=  false, 

objectcnt  :=  1; 

{ 

LOOP  TO  CHECK  EACH  OBJECT  IN  THE  OBJECT  ARRAY. 

} 

repeat 

if  (((picture  objects '[object  cnt|.obj_bsphere  z  >  i  source  *)  and 

ji  ray  i  .»  0  ))  or 

( (picture. objects  '[object  cnt|.obj_bsphere_*  <  i_source_z)  and 
(i  ray  *  <  0  )))  then  begin 

< 

ESTABLISH  PLANE  ON  WHICH  TO  DRAW  BOUNDING  CIRCLE. 

} 

caleq(  picture. objects  '[object _cnt].obj_bsphere_x, 
picture.objects  ‘[object_cnt].obj_bsphere  y, 
piciure.objects  ‘[object  cntj.obj  bsphere  *, 
irayx,  irayy,  i  ray  z, 

A,  B,  C,  D); 

{ 

FIND  INTERSECTION  POINT  BETWEEN  THAT  PLANE  AND  THE  SHOOTING  RAY. 

} 

find  intersection  point ( A, B,C,D, 

i  ray  x,  irayy,  irayz, 
i  source  x,  i  source  y,  i  source  z, 
bplane  intersection  x, 
bplane  intersection  y, 
bplaneintersection  z); 

{ 

DETERMINE  DISTANCE  BETWEEN  THE  CENTER  OF  THE  CIRCLE  AND  THE 
INTERSECTION  POINT 

I 

distance_from_intersection  :  = 

sqrt(sqr(picture. objects  '[object_ent].obj_bsphere_x  - 
bplane  intersection_x)  + 
sqr(picture. objects '[object  cntj.obj  bsphere  y- 
bplane_intersection_y)  + 
sqr(picture. objects '[object  cnt|.obj_bsphere_z  - 
bplane  intersection  z)  ); 


{ 

IF  INTERSECTION  POINT  LIES  WITHIN  CIRCLE  THEN  START  CHECKING  THE 
SUBOBJECTS  THAT  MAKE  UP  THE  OBJECT. 

} 

if  distance  from  intersection  <  = 

picture. objects 'lobject  cntl.obj  bsphere  radius  then  begin 

check_for_subobj  intersection  (object_cnt, 

irayx,  i_ray_y,  i_ray_z, 

i  source  x,  i  source  y,  i  source  z, 

polygonX,  polygonY,  polygonZ, 

subobj  :nt. 

cpart_cnt, 

polygon  cnt, 

found  intersection); 

{ 

DETERMINE  THE  DISTANCE  BETWEEN  THE  ORIGIN  AND  INTERSECTION  POINT  OF 
OF  THE  RAY 
} 

if  found_intersection  then  begin 

view_bplane  vector_x  :=  i_source_x  -  bplane_intersection  x; 
view_bplane_vector_y  :=  i  source  _y  -  bplane  intersection  y; 
view_bplane  vector_z  :=  i_source_z  -  bplane_mtersection  z; 

distancefromviewposition  :=  sqrt(sqr( view  bplane  vector  x)  + 

sqr(view  bplane  vector  y)  + 
sqr(view_bplane_vector_z)); 

{ 

SELECT  ONE  CLOSEST  TO  RAY’S  ORIGIN. 

} 

if  distance_from  viewposition  <  closest  object  then  begin 
closest _object  :=  distance  from  viewposition; 
o_object_idx  :=  object  cnt; 
o_subobj_idx  :=  subobj  cnt; 
o_cpart_idx  :=  cpart  cnt; 
opolyidx  :=  polygoncnt; 
ointersectionx  :=  polygonX; 
o_intersection_y  :=  polygonY; 
o_intersection  z  :=  polygonZ; 
o  intersection  flag  found  intersection 
“nd 
•md 
end; 

end; 

object  cnt  :=  object  cnt  +  1; 
until  (object_cnt  >  picture. num_objs); 

end;  {  check  for  intersection  } 


INTENSITY  PROCEDURES 


{ 


*  . * . . . * . STACK  EMPTY 

*  CALLED  FROM  MAIN  and  POP 

*  CALLS  TO  NONE 


*  DESC  CHECKS  TO  SEE  IF  STACK  EMPTY 

*  INPUT  POINTER  TO  TOP  OF  STACK 

*  OUTPUT  BOOLEAN  VALUE  -  TRUE  FALSE 


function  stackempty  (input _ray_top  :  ray_ptr)  boolean, 
begin 

STACK  EMPTY  :=  input _ray_top  =  nil; 
end;  {  STACK  EMPTY  } 


{ 

. ****** . ***  STACK  EXCEEDED  ****** . * . 

*  CALLED  FROM  :  MAIN 

*  CALLS  TO  :  NONE 

*  DESC  :  CHECKS  TO  SEE  IF  STACK  IS  FULL 

*  INPUT  :  CURRENT  SIZE  OF  STACK  and  MAXIMUM  SIZE  OF  STACK 

*  OUTPUT  :  BOOLEAN  VALUE  TRUE/FALSE 

*************4********************************************************* 

}  function  stack  exceeded(input  tl, 

input_t2  :  integer)  :  boolean; 


begin 

if  inputtl  >  input_t2  then 
stack _exceeded  :=  true 
else 

stack  exceeded  :=  false; 


-n<i:  ■'  STACK  EXCEEDED  } 


{.PA} 

{ 

************************  CALCULATEREFRACTEDRAY  ************************** 

*  CALLED  FROM  :  MAIN 

*  CALLS  TO  :  NONE 

*  DESC  :  This  calculates  the  direction  of  a  refracted  ray. 
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relrai  ;*M  ra.  <•  v»  r»a -<1  i'  n -  »  x-  Xt 


.  r  *  1  u r<*  ii  *eira  ;»•!  *ax  rax  »  *  '•  'r  \ 


i  rax  vector  x 


t  rax  vfclor  i 


i  obj  surface  normal  x 
fbi  »iirf»  -e  n  >rT.  a. 
i  obj  surface  norinii  i  real 
i  tbj  ridx  real 

i  global  ridx  real 
x  ar  o  refracted  ra>  x 
o  refracted  ra>  y 
o  refracted  rax  i  rrai 


vr  ■•Mrs.  i  -x .  as 


testKf. 

testKf2 


i\f  reai. 


Kn  real, 

abs  dot  product  real, 
length  ofray  real, 
v  1  x,  v  1  v,  v  1  i  real. 


THIS  PRODUCES  THE  ABSOLUTE  VALUE  OF  THE  DOT  PRODUCT  OF  THE  INCOMING  RAY 
AND  THE  SURFACE  NORMAL  OF  THE  INTERSECTED  SURFACE 


abs  dot  pt  iduct  :=  Abs((i_obj_surface_normal_x  *  i_ray  vector _x) 
(i  obj  surface _normal  y  *  i  ray  vector  y)  - 
(i  obj  surface_norinal  z  *  i  ray  vector  t)); 


if  abs  dot  product  =  0  then  begi 


r h is  :s  just  \  precu.tion 


o  retracted  ray  dag  .  =  latse, 


o  refracted  ray  _x 
o  refracted  ray  _y 
o  refracted  _ray_x 
end 

else  begin 


1 

THIS  PRODUCES  THE  UNIT  NORMAL  VECTOR  IN  THE  DIRECTION  OF  THE  INCOMING 
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Kn  >bi  r i .  1  x  gc-Gal  ndx 


THE  i  UCI  1  \TION  of  THE  FRESNEL  COEFFICIENT  IS  DIVIDED  I  P  THIS  WAY 
INORDER  To  CHECK  I.  \TFR  FOR  \\  IMAGINARY  DENOMINATOR  WHICH  INDICATES 
TOT  XI  INTERN  AL  REFLECTION 

testKf"  (sqr(Kn)  *  «qr(length  of  ray)) 
t ext  Kf.'  sqr(*qrt  (sqr(i  obj  surface  normal  x  -  vl  x)  - 

■  't'j  -urraor  normal  >  *  vl  y)  - 
»qr(i  oDj  surface  normal  l  *  vl  i  |j). 


.f  1 1 rut  Kf  -  test  Kf2  •  -  0  )  then  begin 

,  i  MAO  IN  aRN  DENOMINATOR  -  TOTAL  INTERNAL  REFLECTION  IS  OCCURING  } 
o  refracted  ray  flag  .  -  false, 
o  refracted  ray  x  -  0: 
o  refracted  ray  y •  :=  0, 
o  refracted  ray  *  =  0 
end 

else  begin 

o  refracted  ray  flag  true; 

Kf  -  1  sqrt(testKf  -  tesiKf2); 

o  refracted  ray  x  :=  (Kf  *  (i  obj  surface  normal  x  •+•  vl  x))  - 
i  obj  surface  normal  x; 

o  refracted  ray  _y  .=  (Kf  *  (i  obj  surface  normal  y  ■+■  vl  y))  - 
i  objsurfacenormaly; 

o  refracted  ray  i  :=  (Kf  *  (i  obj  surface  normal  z  +  vl  z ) )  - 
i  obj  surface  normal  z 

end 

end, 

end  1  calculate  refracted  ray  ) 

\ 

i 

. . .  CALCULATE  REFLECTED  RAY  ************************** 

*  CALLED  FROM  MAIN  and  CALCULATE  INTENSITY 

*  CALLS  TO  NONE 

*  DESC  CALCULATES  THE  DIRECTION  OF  A  REFLECTED  RAY 

*  INPUT  :  Direction  of  the  source  ray  and  surface  normal  of  intersected 

*  object. 

*  OUTPUT  Flag  indicating  existance  of  reflected  ray  If  one  exists 
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then  it’s  direction  is  given. 


************************************************************************** 

} 

procedure  calculate  reflected  ray  (i  ray  vector  x, 

i_ray_vector  y, 
i  ray  _vector  z  :  real; 

i  obj  surface  normal  x. 

i _ obj  surface  normal  y, 

i_obj  surface  normal  z  :  real; 
var  o_reflected_ray_x, 
o_reflected_ray_y, 
o  reflected  ray_.i  ;  real; 
var  o_reflected  rayjflag  :  boolean  ); 


var 

abs_dot  product  :  real; 
length  of_ray  :  real; 
vl  x,  vl_y,  vl  i  :  real; 

Rl_x,  Rl_y.  Rl_z  ;  real; 

begin 

{ 

ABSOLUTE  VALUE  OF  THE  DOT  PRODUCT  OF  THE  OBJECTS  SURFACE  NORMAL  AND 
OF  THE  INCOING  LIGHT  RAY. 

} 

abs  dot  product  :=  Abs((i_obj_surface_normal_x  *  i  ray  vector  x)  + 

(iobjsurfacenormaly  *  i_ray_vector_y)  + 

(*_°bj_surface_normal_z  *  i_ray_vector_z)); 

vl  *  :=  i  ray  vector  x  /  absdotproduct; 
vl_y  :=  i  ray  vector  y  /  abs_dot_product; 
vl_z  :=  i_ray_vector_z  /  abs_dct  product; 

{  CALCULATE  UNIT  NORMAL  VECTOR  IN  THE  DIRECTION  OF  THE  INCOMING  RAY.  } 
Rl_x  :=  vl_x  +  (2  *  i  obj  surface  normal  x); 

Rl_y  :=  vl_y  +  (2  *  i_obj_surface  normal  y); 

Rl_z  :=  vl_z  +  (2  *  i_obj_surface  normal  z); 

length_of_ray  :=  (sqrt(sqr(Rl_x)  +  sqr(Rl_y)  +  sqr(Rl_z))); 

f  CALCULATE  REFLECTED  RAY  } 

>  redecied  rav  •<  —  Rl  x  length  of  ray; 

>  reflected  ray  y  -  R.l  v  length  of  rav; 

o  redected  ray  z  ;=  Rl  z  /  length  oi  ray; 

o_reflected_ray  flag  ;=  true; 

end;  {  calculate_reflected_ray  } 
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{ 

************************  CALCULATE  INTENSITY  ********************** 

*  CALLED  FROM  :  12  and  16 

*  CALLS  TO  :  CHECKFORINTERSECTION  and  CALCULATE  REFLECTED  RAY 

*  DESC  :  Calculates  the  intensity  at  any  given  point. 

*  INPUT  :  Color  component  being  calculated,  current  ray  data, 

*  pointer  to  object  data,  and  pointer  to  light  data. 

*  OUTPUT  :  Intensity  at  a  given  point,  either  to  be  displayed  or 

*  set  in  the  appropriate  source  ray. 

******************************************************************* 

} 

procedure  CALCULATE_INTENSITY  f  input_color  :  colortype; 
input  _ray_d  :  real; 

input_ambient  :  real; 
input_ray_I_t  :  real; 

input_ray_I_s  :  real; 

input  _ray_vector_x, 
input  ray  vector  y, 
input  rav  veccor  i  :  real; 
uiput_number_of_lights  ;  integer; 
input_obj_K  a  :  real; 
input  _obj_K_s  ;  real; 
input  _obj_K  t  ;  real; 
input  obj  K  d  :  real; 
input  _obj_phong  exp  :  integer; 
input  intersection  x  :  real; 
input  intersection  y  :  real; 
input  intersection  s  :  real; 
input_obj_surface_normal_x  :  real; 
input_obj_surface_normal_y  :  real; 
input  obj  surface  normal  z  :  real; 
input  picture  object  :  object  ptr; 
input  light  top  :  light  ptr; 
var  io_intensity  ;  real  ); 

var 

ans  :  char; 

I  d  :  real; 

I  I  :real; 

J,  i  :  integer; 
distance  :  real; 
sight  ray  x, 

-ignt  rav  y. 

^igni  "av  r.  -rai 

unit_sight_x, 

unitsighty, 

unit_sight_z  :  real; 

lightrayx, 

light  ray  _y, 

light  ray  _i  ;  real; 

unitlightx, 

unitlighty, 
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unit_light_z  :  real; 

unit_reflected_x, 

unit  reflected  _y , 

unit_reflected_z  :  real; 

intersection_x, 

intersection  y, 

intersection  z  :  real; 

reflected_Iight_ray_x, 

reflected_light  ray  _y , 

reflected_light_ray_z  :  real; 

source_x, 

source  y, 

sourcez  :  real; 

checkx, 

check_y, 

check  z  :  real; 

obj  idx,  subobjidx,  cpart  idx,  polygon_idx  :  integer; 

intersection  flag  :  boolean; 

obj  light  distance  :  real; 

reflectedray  :  boolean; 

intersected_obj_Kt  :  real; 

begin 

{  THIS  ELIMINATES  THE  SAME  INTERSECTION  POINT  FROM  BEING  SELECTED  AGAIN.} 

if  input  ray  d  >  0.1  then  begin 

I_1  :=~ 0.0; 

iointensity  :=  0.0; 
reflected  ray  :=  false; 
intersectionflag  :=  false; 

{  THIS  SETS  UP  THE  SIGHT  RAY  } 
sightrayx  :=  -input  ray  vector  x; 
sight  ray  y  :=  -input_ray_vector_y; 
sight  ray  z  :=  -input  ray  vector  z; 
distance  :=  sqrt(sqr(sight_ray_x)  + 
sqr(sight_ray_y)  + 
sqr(sight_ray_z)); 

unit_sight_x  :=  sight  ray_x  /  distance; 
unit  sight  v  :=  sight  ray_v  /  distance; 
unit  sight  t,  sight  rav_;  ,  distance; 

for  l  :=  1  to  (input  number_of_lights  +  1)  do  begin 

{  THIS  GENERATES  THE  SHADOW  FEELERS  } 

light_ray_x  :=  picture. lights “ (i) .light _x  -  input_intersection_x; 
light_ray_y  :=  picture. lights “[i].light_y  -  input_intersection_y; 
light  ray  z  :=  picture,  lights  *[i).  light  z  -  input  intersection  z; 
distance  :=  (sqrt(sqr(light_ray_x)  + 
sqr(light_ray_y)  + 


sqr(Iight_ray  _z))); 

{  CONVERTS  IT  TO  A  UNIT  VECTOR  } 
unit_light_x  :=  light_ray_x  /  distance; 
unit_light_y  :=  light_ray_y  /  distance; 
unit  light  _z  :=  light  ray  _z  /  distance: 

source_x  :=  picture. lights  “!i|.light_x; 
source_y  :=  picture. lights  '[ij.light_y; 
sourcez  :=  picture.lights  “[i).light_z; 

{ 

CHECK  TO  SEE  IF  ANY  OF  THE  SHADOW  FEELERS  INTERSECT  ANYTHING. 

} 

check_for_intersection(-unit_light_x,  -unit  light_y,  -unit_light_z, 
sourcex, 
source_y, 
source  z, 

input  pictureobject. 
intersection_x, 
intersectiony, 
intersection  _z, 

objidx,  subobj  idx,  cpart  idx,  polygon  idx, 
in  tersectaon  _ilag) ; 

{ 

CHECK  i O  INSURE  THAT  THE  SAME  POINT  IS  NOT  CONSIDERED  AGAIN,  WHICH  CAN 
HAPPEN. 

} 

checkx  :=  intersection  _x  -  input  intersection  _x; 
check_y  :=  intersection  y  -  input  intersection  y; 
check_z  :=  intersection_z  -  input  intersection  z; 

{ 

PULL  THE  PROPER  CHARACTERISTICS  OF  THE  OBJECT  OUT  TO  DEAL  WITH  THE 
APPROPRIATE  COMPONENT  OF  LIGHT  THAT  IS  BEING  CURRENTLY  DEALT  WITH 

} 

if  color  =  red  then 
I_1  :=  picture.lights‘(i].I_r; 
intersected _obj  Kt  :=  picture. 

objects  ”[obj_idxj. 

sub  objects  “'subobj  idx!. 

common  parts  '  cpart  idxi. 

K  _r.r: 

if  color  =  green  then 
I  J  :=  picture.lights  *[i|.I_g; 
intersected _obj_Kt  :=  picture. 

objects  “[obj_idxj. 
sub_objects  *  [subobj  _idx] . 
common  _parts  *(cpart_idx|. 

K_tg; 

if  color  =  blue  then 


I_1  :=  picture,  lights  *[i].I_b; 
intersected  obj  Kt  :=  picture. 

objects  “[obj  _idx] . 
sub_objects  *[subobj_idxj. 
common_parts  '(cpart_idxj. 

K  _tb: 

obj_light_distance  :=  (sqrt(sqr(intersection_x  - 

picture.lights'[ij.light_x)  + 
sqr(intersection  y  - 

picture.lights'jij. light  y)  + 
sqr(intersection  z  - 

picture,  lights  “(ij.light _ z))); 


{ 

IF  THERE  HAS  BEEN  AN  INTERSETION  AND  THE  POINT  BEING  CONSIDERED  IS  NOT 
THE  ORIGIN  OF  THE  RAY  THEN  CHECK  TO  SEE  IF  THE  OBJECT  IS  OPAQUE.  IF  IT  IS 
THEN  RETURN  TO  BEGINNING  OF  LOOP  AND  CHECK  NEXT  SHADOW  FEELER.  IF  IT  IS 
NOT  OPAQUE  THEN  CALCULATE  INTENSITY  AT  THAT  POINT  AND  CONTINUE  FOLLOWING 
THE  RAY  TO  SEE  IF  IT  INTERSECTS  ANYTHING  ELSE.  CONTINUE  THIS  LOOP 
EITHER  UNTIL  NO  MORE  OBJECTS  ARE  LEFT  OR  UNTIL  AN  OPAQUE  ONE  IS 
INTERSECTED. 

} 

{ 

THIS  PART  OF  THE  CODE  HAS  NEVER  BEEN  TESTED.  I  JUST  CODED  IT  AS  I 
THOUGHT  IT  SHOULD  BE  FROM  THE  ALGORITHM  IN  ROGERS  BOOK  PP.  377. 

QUITE  FRANKLY  I  STILL  DON’T  FULLY  UNDERSTAND  WHAT  IS  SUPPOSE  TO  TAKE  PLACE 
HERE. 

} 

if  intersection  flag  then  begin 
while  ((intersectionflag)  and 
((check  x  >  0)  or 
(check_y  >  0)  or 
(check  x  >  0)))  do  begin 


intersection_flag  :=  false; 


if  not(  intersected  obj  Kt  =  0)  then  begin 
if  (input_color  =  red)  then  begin 

I_1  :=  picture. lights  “[i].I_r  * 
intersected  obi  Kt: 

f  input  :olor  =  ^r^eni  fien  besjm 

I_1  :=  picture. lights "(i].I_g  * 
intersectedobjKt; 

end; 

if  (input  color  =  blue)  then  begin 

I_1  :=  picture. lights *[i].I_b  * 
intersected_obj  Kt; 
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end; 

calculate _reflected  ray  (input  _ray  _vector  _x, 
input  _ray  _vector_y , 
input_ray_vector_z, 
input  _obj  surface  jiormal  x, 
input  _obj  surtace_norTnai_y, 
input  _obj_surface_normal_z, 
reflected  ray  x, 
reflected  ray_y, 
reflected  ray_z, 
reflected  ray); 

distance  :=  sqrt(sqr(reflected_ray_x)  ■+■ 
sqr(reflected_ray_y)  + 
sqr (reflected  ray  _z ) ) ; 

unit  reflec ted _x  :=  reflected  ray _x  /  distance; 
unit  reflected  y  :=  reflected  ray_y  I  distance: 
unit  reflected_z  :=  reflected _ray_z  /  distance; 


io  intensity  :=  io  intensity 

((I_l  *  input_obj_K_d)  * 

((inputobjsurfacenormalx  *  unitlightx)  + 
(input_obj_surface_normal_y  *  unitlighty)  + 
(input  obj  surface  normal  z  *  unitlightz)))  + 
((I  I  *  input  obj_K_s)  * 

((unit  sight  x  *  unit  reflected _x)  + 

(unit  sight  y  *  unitreflectedy)  + 

(unitsightz  *  unitreflectedz))); 

source_x  :=  intersection_x; 
source  y  :=  intersection_y; 
source_z  :=  intersection_z; 

check_for_intersection(-unit_light_x, 

-unit_light_y, 

-unitlightz, 
source_x, 
source_y, 
source  z. 

input  picture  object, 
intersection  x, 

intersection_y, 
intersection  z, 
obj  idx,  subobj_idx, 
cpart_idx,  polygon  idx, 
intersectionflag); 

check  x  :=  intersection_x  -  input  intersection  x; 
check  y  :=  intersection  _y  -  input -intersection  y; 
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check  *  :=  intersection _z  -  input  intersection  *; 

end  {if} 
end  {while} 
end  {if} 
else  begin 

calculate  _rerlected  -ay( input  ray  vector  x. 

input  _ray_vector_y, 

input_ray_vector_z, 

input  _obj_surface_normal_x, 

input_obj_surface_noriral_y, 

input  _obj  surface  normal  *, 

reflected  _ray_x, 

reflected  _ray_y, 

reflected_ray_*, 

reflected  ray); 

distance  :=  sqrt(sqr(reflected_ray_x)  + 
sqrfreflected  ray  y) 
sqr(retiec  ted  _ray  _z ) ) ; 

unit  reflected _x  :=  reflected  ray _x  /  distance; 
unit  reflected_y  :=  reflected  _ray_y  /  distance; 
unit_reflected_z  :=  redected_ray_z  /  distance; 

iointensity  :=  iointensity  + 

((I_l  *  input_cbj_K_d)  * 

((inputobjsurfacenormalx  *  unitlightx)  + 

(inputobjsurfacenormaiy  *  unitlight  y)  + 

(inputobjsurfacenormalz  *  unit  light  *)))  + 

{(I  I  *  inputobjKs)  * 

((unitsightx  *  unit  reflected _x)  + 

(unit  sight  y  *  unitreflectedy)  + 

(unit  sight  *  *  unit  reflected  z))); 
end  {  **  ELSE  **} 
end;  {  **  FOR  **  } 

{ 

THIS  IS  THE  STUB  TO  JUST  HAVE  EVERY  OBJECT  ILLUMINATED  BY  AMBIENT  LIGHT 
io  intensity  :=  input  obj  K  a  *  input  ambient; 


f  I******************************************************  | 

THIS  iS  WHERE  THE  FINAL  INTENSITY  IS  CALCULATED 

i 

io_intensity  :=  ((input  obj_K  a  *  input_ambient)  + 
iointensity  + 

((input  obj  K  s  *  input  ray _I_s)  + 

(input_obj_K_t  *  input _ray_I_t)))  ; 

{/inputray  d  or  /2  or  / 1 } 


( 


THIS  IS  JUST  TO  KEEP  ALL  VALUES  WITHIN  A  RANGE  WHERE  THEY  CAN  BE  DISPLAYED 

if  io  intensity  >1.00  then 
io_intensity  :  =  1.00; 
if  io  intensity  <  0.00  then 
io  intensity  :=  0.00; 


{  .*•*».«.*•«**********•******,*«***»««»•».«,***««***«****«  | 
end 

else  begin 

{ 

IF  THE  INPUT  RAY  DISTANCE  IS  LESS  THEN  ONE  THAN  YOU  ARE  CONSIDERING  THE 
SAME  POINT  AND  HENCE  THE  INTENSITY  THERE  SHOULD  BE  0. 

} 

io  intensity  :=  0.0 
end; 


end;  {  CALCULATE  INTENSITY  } 

{.PA} 

{ 

. . DISPLAY  PIXEL  *******.•*****.***•******»*»* 

*  CALLED  FROM  :  MAIN 

*  CALLS  TO  :  NONE 

*  DESC  :  WRITES  OUTPUT  TO  FILE 

*  INPUT  :  THE  VALUES  FOR  THE  RED,  GREEN,  AND  BLUE  COMPONENTS  OF  LIGHT. 

*  OUTPUT  :  NONE 

It******************************************************************* 

} 

procedure  DISPLAY  PIXEL  (  inputintensityred, 
input  intensity  green , 
input  intensity  blue  :  real; 
input  jpixel_x, 
inputpixely, 
input_pixel_z  :  real  ); 


begin 

write  (outfile, input  _intensity_red:3:2); 
write  (outfile,’  ’,input_intensity_green:3:2); 
writeln  (outfile,’  ’,input_intensity_blue:3:2); 

■nu:  DISPLAY -PIXEL  } 

V.PA} 

{ 

***** . . . **********  POP  ******* . * . * . ***** 

*  CALLED  FROM  :  MAIN,  12,  16 

*  CALLS  TO  :  NONE 

*  DESC  :  Removes  a  ray  from  the  top  of  the  stack. 

*  INPUT  :  Pointer  to  the  current  top  top  of  stack. 

*  OUTPUT  :  The  ray  just  popped  from  the  stack  and  a  pointer  to  the  new  top 
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*  top  of  stack. 

************************************************************************* 

} 

procedure  POP  (var  outputraytype  :  raytype; 

var  output_ray_origin_x  :  real; 

var  output  _ray  origin  y  ■  real; 

var  outpuc_ray_origin  z  :  real; 

var  output_ray_vector_x  :  real; 

var  output_ray_vector_y  :  real; 

var  output_ray_vector  z  ;  real; 
var  output _ray_stype  :  raytype; 

var  output  intersection  flag  .  boolean; 

var  outputobjidx  :  integer; 

var  output  subobj  idx  :  integer; 
var  output_cpart_idx  ;  integer; 

var  outputpolygonidx  :  integer; 

var  output_intersection  x  ;  real; 

var  output  intersection  _y  :  real; 

var  output  intersection  r,  :  real; 
var  outputd  :  real; 

var  output_I_tr,  output_I  tg,  output  I  tb  ;  real; 
var  output_I_sr,  output_I  sg,  output  _I_sb  :  real; 
var  io  top  :  rayptr  ); 

begin 


if  (stack  empty  (iotop))  then  begin 
writeln(’STACK  UNDERFLOW  ERROR’) 
end  {  if  } 
else  begin 

output  ray  type 
output  ray  origin  x 
output  ray  origin  y 
output  ray_origin_z 
output_ray_vector_x 
output  ray  vector_y 
output  ray_vector_z 
output  ray  stype 


:=  iotop  .ray  type; 

;=  iotop  ".rayoriginx; 
:=  io  top ‘.ray  origin_y; 
;=  io_top  *.ray_origin_z; 
;=  io_top '.ray  vector  x; 
:=  io_top ‘.ray  vector  y; 
;=  io_top‘.ray  vector_z; 
:=  io  top ‘.ray  stype; 


output  intersection  flag  :=  io  top ‘.intersection  flag; 


output_obj  idx 
output  subobj  _idx 
>utt>uc  'parr,  idx 
nitpuc  polygon  idx 
■jiupuc  intersection  x 
output  intersection  y 
output  intersection_z 
output  d 
output  I  tr 
output  1  tg 
output_I_tb 
output  I  sr 
output  I  sg 


:=  io_top ‘.obj_idx; 

:=  io  top  ‘.subobj  idx; 

—  io  _fop  *  ;oart  idx; 

•  =  so  op  poivgon  dx; 
—  ;o  ion  '  .n  terser  non  x: 
:=  io_top ‘.intersection  y; 
:=  io_top ‘.intersection  z; 
:=  io_top  ".d; 

:=  io  top‘.I_tr; 

:=  io_top  *.I_tg; 

:=  io_top‘.I_tb; 

;=  io  top ‘.I_sr; 

:=  io_top  ‘.I_sg; 
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* 


'In 


output_I_sb  :=  io  top  M  sb; 

r&y_next  :=  io_top *. ray  link; 

{  remove  old  pointer  } 
dispose(io_top); 

{  set  top  of  stack  pointer  to  new  top  of  stack  } 
io_top:=ray  next; 
end;  {  eise  } 
end;  {  POP  } 


{PA} 


**************************** 


PUSH 


****************** 


*  CALLED  FROM:  MAIN,  12,  16 

*  CALLS  TO  :  NONE 

*  DESC  :  Places  a  ray  on  the  top  of  the  stack. 

*  INPUT:  The  current  top  of  the  stack,  and  the  data  for  a  new  ray 

*  OUTPUT:  The  pointer  to  th  new  top  of  stack. 


procedure  PUSH  (  input  ray  type  :  ra) 

input_ray_origin  x  :  real; 
input  ray _origin  y  :  real; 
input  ray  origin  z  :  real; 
input_ray  vector  x  :  real; 

input  ray  vector  y  :  real; 

inputrayvectorz  :  real; 
input  ray  stype  :  raytype; 

input  intersection  flag  :  boolean; 
input  obj  idx  ;  integer; 

inputsubobjidx  :  integer; 
inputcpartidx  :  integer; 

inputpolygonidx  :  integer; 

input_intersection  x  :  real; 
inputintersection  y  :  real; 
input  intersection  z  :  real; 


raytype; 


input  obj  idx  :  integer, 

input  subobj  idx  :  integer; 
input  cpart  idx  :  integer; 
input  polygon  idx  :  integer; 
input_intersection_x  :  real; 
input  intersection  y  :  real; 
input_intersection  z  :  real; 
inputjd  :  real; 

input  l  tr,  input  l  tg,  input_I_tb  :  real; 
input  l  sr,  input  l  sg,  input_I_sb  :  real; 
var  io_top  :  ray  ptr  ); 


newiray  rtirrent,): 
ray  --urrent  ray  type 

rayjcurrent  .ray_origin_x 
ray_current  ".ray  origin  y 
ray_current  ‘.ray  origin  z 
ray_current  \ray_vector_x 
ray_current  *.ray_vector_y 
ray_current ‘.ray  vector  z 
ray_current ‘.ray  stype  := 
ray_current  “,intersection_flag  : 


input  ray  type: 


:=  input_ray  origin_x; 

:=  input_ray_origin_y; 

=  input_ray_origin_z; 

:=  input  ray  vector_x; 

:=  input  ray  vector_y; 

:=  input_ray  vector_z; 

=  input  ray  stype; 

:=  input  intersection  flag; 


ray  current 
ray  current 
ray  current 
ray_current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
ray  current 
raycurrent 
io  top 


obj_idx  :=  input  obj_idx; 

subobj_idx  :=  input  jsubobj  idx; 
.cpartidx  :=  input  cpart  idx; 

.polygon_idx  :=  input_polygon_idx; 

. intersection _x  :=  input  intersection  x; 

intersection  v  —  input  intersection  v; 

intersection^  input  intersection  v, 

.d  :=  input  d; 

.I_tr  :=  input_I_tr; 

•  I_tg  :=  input_I_tg; 

.I_tb  :=  input  I  tb; 

I_sr  =  input  I  sr; 

•  I_sg  :=  input  _I_sg; 

.Ijsb  :=  inpuc_I_sb; 

'.ray_link  :=  io_top; 

:=  ray  current; 


end;  {  PUSH  } 


MAIN 


I 


A 

i 


$ 

i 


i  hT 


PROGRAM  RAYTRACER; 

{ 

44444444444444444444444444444444444444*44*44444444*4444444444*444444444 

«*-«*******'«*K***ac**»:**»*«******a:**»»**K**^-»:»***x**jt**x******»****)K*«*** 

**  PROG  :  RAY. PAS 
**  AUTHOR  :  Paul  G.  Smith 
**  DATE  :  11  May  1987 

**  DESC  :  A  ray  tracing  prototype  with  a  global  illumination  model 
'*  integrated  into  it. 

**  INPUT  :  A  sequential  scene  file  under  PICTURE. PAS. 

**  OUTPUT  :  A  bitmap  file  containing  the  red,  green,  and  blue  color 
**  color  components.  Their  values  range  from  0-1  and  need  to 

**  be  converted  for  display  on  an  RGB  color  monitor. 

****************************************************************4****** 

44*44*44*44*4444444***4**44444*44*444**4444444444**44444444444444444444 


{  *  INCLUDE  FILES  *  } 

{SI  declared. pas!  -  {  DECLARATION  SECTION  } 

{.PA} 

{SI  intprcsft.pas}  -  {  INTERSECTION  PROCEDURES  } 

{.PA} 

{$1  procsd.pas}  _  {  INTENSITY  and  UTILITY  PROCEDURES  } 


{.PA} 

{ 

***********************************  j2  ******************************** 


*  CALLED  FROM:  MAIN 

*  CALLS  TO:  PUSH,  POP,  CALCULATE  INTENSITY 

*  DESC  :  Calculates  the  light  intensity  at  a  given  intersection  point 

*  INPUT  :  A  complete  ray  data  record,  the  pointer  to  the  light  array 

*  and  the  pointer  to  the  object  array 

*  OUTPUT  :  The  intensity  at  a  given  intersection  point.  If  the  input 

4  was  a  view  ray  then  this  intensity  will  be  the  intensity 

'  displayed.  If  the  input  ray  is  a  reflected  ray  then  this 

*  intensity  is  assigned  to  the  I  s  field  in  the  source  ray. 

*  If  the  input  ray  is  a  refracted  ray  then  this  intensity  is 
assigned  fo  the  I  t  field  in  the  source  rav. 


procedure  12  (i  d  :  real; 

i _ I _ tr,  i_I_tg,  i_I_tb  :  real; 

i_I^sr,  i_I_sg,  i_I_sb  :  real; 
i_number  of  light  sources  :  integer; 
i_ambient  r,  i  ambient  g,  i  ambient  b  :  real; 
i_K  ar,  i_K  ag,  i  K_ab  :  real; 
i_K_sr,  i_K  sg,  i_K_sb  :  real; 
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real; 

:  real; 


real; 


i_K_tr,  i_K_tg,  i_K_tb 
i_K_dr,  i_K_dg,  i_K_db 
i_ray_type  ;  raytype; 
i_ray_top  :  ray_ptr; 
i_light_top  :  light  ptr; 
i  _obj  _ptr  :  object  ptr; 
i_obj_pnong_exp  ;  integer: 
i_ray_vector_x, 
i_ray_vector_y, 
irayvectorz 
i_intersection_x, 
i_intersection  y. 
i_intersection  z 
i_surface_normal_x, 
i_surface  normal  y, 
i_surface_normal_z 
var  io_ray_generation  number  :  integer; 
var  o_intensity_red, 
ointensity  green, 
o  intensity  blue  :  real  ); 


real; 


;  real; 


{  SET  UP  TEMPORARY  AREA  TO  HOLD  RAYS  POPPED  FROM  STACK  } 

templ_ray  type  ;  raytype; 


real; 


:  real; 
raytype; 


tempi  ray  origin  x, 
tempi  ray  originy, 
templrayoriginz  : 
tempi  ray  vector  x, 
tempi  ray  vector  y, 
tempi  ray  vector  z  : 
tempi  ray  stype 
templ_intersection_flag  ;  boolean; 
templ_obj  idx, 
templsubobjidx, 
temp l  cpart  idx, 
templ_polygon_idx 
tempi  intersection  x, 
tempi  -intersection  y, 
templintersectionz 
templ  d  :  real; 
tempi  I_tr,  tempi  I_tg,  templ  l 
temDl  I  sr.  tempi  l  sg,  tempi  I 


:  integer; 


real; 


tb 

\b 


real; 

real: 


remo2  ray  'ype  raytype: 

tempi  ray  origin  x, 
temp2_ray_origin_y, 
temp2_ray_origin  z  :  real; 
temp2_ray  vector  x, 
temp2_ray  vector  y, 
temp2_ray _vector_z  :  real; 
temp2 _ray_stype  :  raytype; 
temp2_intersection  flag  :  boolean; 


"*  *  -w  r . n 't. .a  Tj  r ia 


temp2  obj  idx, 
temp2  subobj  idx, 
temp2_cpart  idx, 
temp2_polygon  idx  :  integer; 
temp2_intersection_x, 
temp2  intersection  v. 
cemp2  intersection  z  :  real; 
temp2_d  :  real; 

temp2_I_tr,  temp2_I_tg,  temp2_I_tb  :  real; 
temp2_I_9r,  temp2 __I_sg,  temp2_I_sb  :  real; 

begin 

{  CALCULATE  INTENSITY  OF  THE  RED  COMPONENT  OF  LIGHT  } 

calculate_intensity(  red, 

Ld- 

i  ambient  r, 
i  I  tr, 
i  I  sr. 

i_ray_vector_x,  i_ray  vector_y,  i  ray  vector  z, 

i_number_of_light  sources, 

i_K_ar, 

i_K_sr, 

i  K  tr, 

i_K_dr, 

i_obj  phong  exp, 

iintersectionx, 

iintersectiony, 

iintersectionz, 

isurfacenormalx, 

i_surface_normal_y, 

isurfacenormalz, 

iobjptr, 

i  Jighttop, 

o  intensity_red  ); 

{  CALCULATE  INTENSITY  OF  THE  GREEN  COMPONENT  OF  LIGHT  } 
calculate_in tensity (  green, 
id, 

i_ambient_g, 
i_I_tg, 
i  I  sg, 

;  rav  vector  x,  i  ray  vector  v,  i  rav  vector 

;  number  of  light  sources. 

i_K_ag, 

i_K_sg, 

i_K_tg, 

i_K_dg, 

i  obj_phong_exp, 
i  intersection  x, 
i  intersection  y, 
i  intersection  z, 
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i_surface  normal  x. 
i_surface_normal  y, 
i  surface  normal  *. 
i_obj_ptr, 
i  light  top, 
o  intensity  green  1; 

{  CALCULATE  INTENSITY  OF  THE  BLUE  COMPONENT  OF  LIGHT  } 
calculate_intensity(  blue, 

i_d, 

i_ambient_b, 

i_I_tb, 

i_I_sb, 

i_ray_vector  x,  i  ray  vector  y,  i_ray_vector_z, 
i_numb«r_of  light_sources, 
i_K_ab, 
i  K  sb, 
i  K  tb, 

i  K  jib. 

lobjphongexp, 
i  intersection  x, 
i  intersection  y, 
i  intersection  z. 

1  surtace_normal  x, 
i  surface  normal  y, 
i  surface  normal  z, 
iobjptr, 
i  light  top, 
o  intensity  blue  ); 

if  (i  ray  type  =  view)  then  begin 

{ 

THE  VIEW  RAY  IS  ALWAYS  THE  LAST  RAY  ON  THE  STACK  WHEN  IT  IS  POPPED  THE 
INTENSITY  DETERMINED  FOR  IT  IS  PASSED  BACK  INTO  MAIN  FOR  DISPLAY. 

} 

{  nothing  } 
end 

else  begin 

if  (i  ray  type  =  reflected)  then  begin 

r 

SINCE  THIS  iS  THE  REFLECTED  HaY  THEN  JUST  ONE  RAY  NEEDS  TO  BE  POPPED  TO 

<;  \;n  tor  ess  to  the  source  ray 

} 

pop  (tempi _ray_type, 

tempi  ray _origin_x,  tempi  ray _origin_y,  tempi  ray _origm_z, 

templrayvectorx,  tempi  ray  vector  y,  tern  pi  ray  vector*, 

tempi  ray  stype, 

tempi  intersection  flag, 

tempi  obj  idx, 

tempi  subobj  ldx, 
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tempi  cpart  idx, 
tempi  polygon  idx, 

tempi  intersection _x,  tempi  _intersection_y .  tempi  intersection  z, 
temp  I  d, 

templ_I_tr,  templ_I_tg,  templltb, 
tempi  I_sr.  tempi  I_sg,  templ_I  sb. 

ray  top  );  - 

{  SET  INTENSITY  IN  SOURCE  RAY  } 
tempi  I  sr  :=  o  intensity  red; 
tempi  I_sg  :=  o  intensity  green; 

tempi  I  sb  :=  o  intensity  blue; 

{  RESTORE  STACK  } 

push(  tempi  ray  type, 

tempi  ray  origin  x,  templ_ray_origin_y,  tempi  ray  origin  z, 

templrayvectorx,  tempi  _ray_vector_y,  tempi  _ray_vector_z, 

tempi  ray  stype, 

tempt  intersection  Jlag, 

templ_obj  idx, 

templ_subobj_idx, 

tempi  cpart  idx, 

tempi  polygon  idx, 

tempi  intersection_x,  tempi  intersection  y,  tempi  intersection  z, 
templd, 

templ  l  tr,  templ  l  tg,  templ  l  tb, 
templ  l  sr,  templ  l  sg,  templ  l  sb, 
ray  top  ); 


else  begin 

{ 

SINCE  THIS  IS  THE  REFRACTED  RAY  TWO  RAYS  MUST  BE  POPPED  TO  GAIN  ACCESS 
TO  THE  SOURCE  RAY. 

} 

pop  (tempi  ray  type, 

templ  ray  origin  x,  temp l_ray  origin _y,  templ  ray  origin  z, 

templ_ray  vector_x,  tempi  ray  vector  y,  tempi  ray  vector  z, 

tempi  ray  stype, 

tempi  intersection  flag, 

templ_obj  idx, 

tempi  subobj  idx. 

rempl  '■part  idx, 

temp l  polygon  idx. 

templ_intersection_x,  tempi_mtersection  y,  tempi  intersection  z, 
templd, 

templ_I_tr,  templ_I_tg,  templ  l  tb, 
tempi  I_sr,  templ_I_sg,  templ_I_sb, 
ray_top  ); 

pop  (temp2_ray_type, 

temp2_ray_origin_x,  temp2_ray_origin_y,  temp2_ray  origin _z, 
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temp2_ray_vector_x,  temp2_ray  ^vector  y.  temp2  ray  vector  z, 

temp2_ray  stype, 
temp2_intersection_flag, 
temp2_obj_idx, 
temp2  jsubobj  _idx, 
temp2  cpart  idx, 
temp2_poiygon_idx, 

temp2_intersection_x,  temp2_intersection_y,  temp2_intersection_z, 
temp2_d, 

temp2  I_tr,  temp2_I_tg,  temp2_I_tb, 
temp2_I_sr,  temp2_I_sg,  temp2_I_sb, 
ray_top  ); 

{  SET  INTENSITY  IN  THE  SOURCE  RAY.  } 
temp2_I_tr  :=  o  intensity  red; 
temp2_I_tg  :=  o  intensity _green; 
temp2_I_tb  :=  o_intensity_biue; 

{  RESTORE  STACK.  } 

push(  temp2_ray_type, 

temp2_ray_origin_x,  temp2_ray_origin_y,  temp2_ray_origin_z, 

temp2  rayvectorx,  temp2_ray_vector_y,  temp2  ray  vectorz, 

temp2  ray  stype, 

temp2  intersection  flag, 

temp2  obj  idx, 

temp2  subobj  idx, 

temp2_cpart_idx, 

temp2_polygon_idx, 

temp2_intersection_x,  temp2_intersection_y,  temp2_intersection_z, 
temp2_d, 

temp2_I_tr,  temp2_I_tg,  temp2_I_tb, 
temp2_I_sr,  temp2_I_sg,  temp2_I_sb, 
ray_top  ); 

push(  templ  ray  type, 

templrayoriginx,  templ  ray  origin  y,  tempi  ray  origin  z, 

templ_ray_vector_x,  tempi  ray _vector_y,  templ  ray  vector  z, 

tempi  ray  stype, 

tempi  intersection  flag, 

templobjidx, 

tempi  subobj  idx, 

tempi  cpart  idx. 

rempl  polygon  idx, 

rempl  intersection  x.  rempl  intersection  v,  tempi  intersection 

templd, 

templ_I_tr,  tempi  I_tg,  templ_I_tb, 
templlsr,  templ_I_sg,  templ_I_sb, 
ray  _top  ) ; 

end; 

{ 

SINCE  THE  RAY  SENT  INTO  THIS  PROCEDURE  IS  NO  LONGER  NEEDED  IT  IS 


DISCARDED,  HENCE  THE  RAY  COUNT  NEEDS  TO  BE  DECREMENTED. 

} 

io  ray  generation  number  :=  io_ray  generation  _number  -  1; 
end; 

end; 

{.PA} 

{ 

******************.*******.*«****  I6  .,*.,*«****•*,***,,***«,,,.***,**, 

*  CALLED  FROM  :  MAIN 

*  CALLS  TO  :  PUSH,  POP,  CALCULATE  INTENSITY 

*  DESC  :  Used  to  calculate  intensity  at  node  of  stack  storage  is 

*  exceeded. 

*  INPUT  :  Incomplete  ray  data  record,  this  ray  can  not  be  continued 

*  because  there  is  no  room  on  the  stack  for  it.  Also  pointers 

*  to  the  light  array  and  object  array. 

*  OUTPUT  :  Intensity,  I  t  or  I_s,  is  set  in  source  ray. 

********************* *^ ************************************************ 

} 

procedure  16  (i  d  :  real; 

i_I_tr,  i_I_tg,  i_I_tb  :  real; 
i_I_sr,  i_I_sg,  i_I_sb  :  real: 
i_number_of_light_sources  :  integer; 
i  ambient  r,  i  ambient  g,  i  ambient  b  :  real; 
i_K_ar,  i_K_ag,  i_K_ab  :  real; 
i_K_sr,  i_K_sg,  i_K  sb  :  real; 

i_K_tr,  i_K_tg,  i_K_tb  :  real; 

i~K_dr,  i_K_dg,  i_K_db  :  real; 
i  ray  type  :  raytype; 
i  ray  top  :  ray  ptr; 
i  light  top  :  light  ptr; 
i  obj  ptr  :  object  _ptr; 
i  obj  phong  exp  :  integer; 
i  ray_vector_x, 
i  ray_vector_y, 
i  ray  vector  z  :  real; 
iintersectionx, 
i  intersection_y, 
i  intersection^  :  real; 
i_surface_normal_x, 
i  surface  normal  y, 
i  ^urtace  normal  c  real; 
var  o  intensity  red. 
o  intensity  green, 
o  intensity  blue  :  real  ); 


var 

{  STORAGE  FOR  POPPED  RAYS  } 
tempi  ray  type  :  raytype; 
tempi  _ray_origin  x, 
tempi  ray _origin_y, 
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tempi  ray  origin_z  :  real: 
templ_ray_vector_x, 
tempi  ray  vector_y, 
templ_ray_vector_z  :  real; 
templ_ray_stype  :  raytype; 
tempi  intersection  flag  ■  boolean: 
tempi  obj  idx. 
tempi  subobj  idx, 
tempi  cpart  idx, 
templ_polygon_idx  :  integer; 
temp  1  intersection  _x, 
tempi  intersection  y, 
tempi  intersection  s  :  real; 
tempi _d  :  real; 

tempi  I  tr,  tempi  I_tg,  templltb  :  real; 
templ  l  sr,  tempi  I_sg,  templ_l_sb  :  real; 

temp2  ray  type  :  raytype; 
temp2  ray  origin 
temp2_ray_origin_y, 
temp2  ray  origin  z  :  real; 

temp2  ray  vector  x, 

temp2  ray  vector  y, 

temp2_ray_vector_i  :  real; 
temp2_ray_stype  :  raytype; 
temp2_intersection_flag  :  boolean; 
temp2_obj_idx, 
temp2_subobj_idx, 
temp2  cpart  idx, 
temp2_polygon  idx  :  integer; 
temp2  intersection _x, 
temp2_intersection_y, 
temp2_intersection_z  :  real; 
temp2_d  :  real; 

temp2_I_tr,  temp2_I_tg,  temp2_I_tb  :  real; 
temp2_I_sr,  temp2_I_sg,  temp2_I_sb  :  real; 

begin 

{***  COULD  INSERT  TREE  EXTENSION  PROCEDURE  ***} 

f  CALCULATE  RED  COMPONENT  OF  LIGHT  } 

calculate  intensity!  red. 

i  *i. 

i  ambient  r, 

i_I_tr, 

i_I_sr, 

i_ray_vector_x,  i_ray_vector_y,  i  ray  vector  z, 
i_number_of_light  sources, 
i  K  ar, 


I  ljs._dr. 

|  i_obj  phong _exp, 

I  i  intersection _x, 

i_intersection_y, 
i  intersection  z, 

I  i  surface  normal  _x, 

I  i  j>urface_normai_y. 

I  i  surface_normal  z. 

i_obj_ptr, 
i_light_top, 
o_intensity_red  ); 
i 

;  {  CALCULATE  GREEN  COMPONENT  OF  LIGHT  } 

calculate_intensity(  green, 

i  J, 

i_ambient_g, 

ijjg, 

i_I_s8i 

i  ray  vector  x.  i  ray  vector  y,  i  ray  vector  z. 
i  i_number_of  iignt  sources, 

;  i_K_ag, 

I  i_K_sg, 

i_K_dg, 

i  obj  phong  exp, 

iintersectionx, 

iintersectiony, 

i  iintersectionz, 

isurfacenormalx, 

isurfacenormaly, 

isurfacenormalz, 

iobjptr, 

ilighttop, 

ointensity  green) ; 

{  CALCULATE  BLUE  COMPONENT  OF  LIGHT  } 
calculate  intensity(  blue, 
id, 

i  ambient  b, 

ijtb, 

i_I_sb, 

i  ray  vector  x.  i  ray  vector  v.  i  ray  vector  z, 
i  numoer  )i'  light  sources, 
i  K  ib. 

i_K_sb, 

i_K_tb, 

i_K_db, 

i_obj_phong  exp, 
iintersectionx, 
i  intersection  y, 
i  intersection  z, 
i  surface  normal  x, 
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i  surface  normal_y, 
i_surface_normal_z, 
iobjptr, 
i_light_top, 
o_intensity  blue) ; 

if  (i  ray  type  =  reflected)  then  begin 

{ 

IF  INPUT  RAY  IS  A  REFLECTED  RAY  THEN  ONE  RAY  MUST  BE  POPPED  TO  GAIN 
ACCESS  TO  IT’S  SOURCE  RAY. 

} 

pop  (templraytype, 

tempi  ray  origin _x,  templ_ray_origin_y,  tempi  _ray_origin_z, 

tempi jray_vector_x,  temp  1  ray _vector_y,  templ_ray_vector_z, 

tempi  ray  stype, 

tempi  intersection  flag, 

tempi  _obj  _idx, 

tempi  subobj  idx. 

tempi  cpart  idx, 

tempi  _polygon  _idx, 

templintersectionx,  templ  intersection  y,  tempi  intersection  s, 
tempi  _d, 

tempi  I  tr,  templltg,  templ  l  tb, 
templlsr,  templlsg,  templ  l  sb, 
i  ray  top  ); 

{  SET  I  s  IN  SOURCE  RAY  } 
templ  l  sr  :=  o  intensity  red; 
tempi _I_sg  :=  o  intensity  green; 
tempi _I_sb  :=  o  intensity  blue; 

{  RETORE  STACK.  } 

push(  templ  ray  type, 

tempi _ray_origin_x,  templ_ray_origin_y,  tempi  _ray_origin_z, 

tempi  ray_vector_x,  templ_ray_vector  y,  templ_ray  vector_z, 

templ_ray  stype, 

tempi  intersection  flag, 

templobjidx, 

templsubobjidx, 

tempi  cpart  idx, 

tempi  oolvgon  idx. 

:pmol  iniersecuon  tempi  intersection  v,  tempi  intersection 
tempi  ii. 

tempi  I_tr,  templ_I_tg,  templ_l_tb, 
templ  l  sr,  templ  l  sg,  templ  l  sb, 
ray  top  ); 


IF  INPUT  RAY  IS  A  REFRACTED  RAY  THEN  TWO  RAYS  MUST  BE  POPPED  FROM  THE 
STACK  TO  GAIN  ACCESS  TO  IT’S  SOURCE  RAY. 

} 

pop  (tempi  _ray_type, 

templ_ray_origin_x,  templray _origin_y,  templ  ray  origin  z, 

tempi  ray  vector _x.  tempi  ray  vector  y.  tempi  ray  vectu.  z. 

tern  ./j.  ray  _stype, 

tempi  _intersection_flag, 

templobjidx, 

tempi  subobj_idx, 

tempi  _cpart  _idx, 

tempi  polygon  idx. 

tempi  intersection _x,  tempi  intersection  y,  tempi  intersection  z, 
templd, 

templ  l  tr,  templltg,  templltb, 
templlsr,  templlsg,  templlsb, 
ray  top  ); 

pop  (temp2  ray  type. 

temp2_ray_origin_x,  temp2_ray _origin_y,  tern  p2_ray  origin  z, 

temp2_ray_vector_x,  temp2_ray_vector_y,  temp2  ray  vector  z, 

temp2_ray_stype, 

temp2_intersection  Sag, 

temp2_obj_idx, 

temp2_subobj_idx, 

temp2_cpart_idx, 

temp2  polygon  idx, 

temp2_intersection_x,  temp2_intersection  y,  temp2  intersection  z, 
temp2_d, 

temp2_I_tr,  temp2_I_tg,  temp2_I_tb, 
temp2_I_sr,  temp2_I  sg,  temp2_I_sb, 
reytop  ); 

{  SET  THE  I  t  FIELD  IN  THE  SOURCE  RAY  } 
temp2_I  tr  :=  ointensity  red; 
temp2_I  tg  .=  o  intensity  green; 
temp2_I_tb  :=  o_intensity_blue; 

{  RESTORE  THE  STACK  } 
push(  temp2_ray_type, 

temp2_ray_origin_x,  temp2_ray  origin  y,  temp2_ray_origin_z, 

•  emp2  *av  vector  x.  temo2  rav  vector  v  •.emo2  ray  vector  :. 

"mp2  -av  'tvp». 
emo2  ntersecrion  .lair. 

temp2_obj  idx, 
temp2  subobj  idx, 
temp2  cpartidx, 
temp2  polygon  idx, 

temp2  intersection  x,  temp2  intersection  y,  temp2  intersection  z, 
temp2_d, 

temp2_I_tr,  temp2  I_tg  temp2  I  tb. 
temp2  I  sr,  temp2  I  sg,  temp2  I  sb, 


ray  sop 


push(  tempi  _ray_type. 

templ_ray_origin_x,  templ_ray_origin_y,  tempi_ray_origin_z, 

templ_ray  origin  x,  templ_ray  origin  y,  tempi  ray  origin  z, 

tempi  ray  stvpe. 

temp  l  _interseccion  riag, 

templ_obj_idx, 

temp  1  _subobj  _idx, 

templcpartidx, 

tempi  polygon  idx, 

tempi  intersection  x,  temp  ^intersection  y,  tempi  intersection  z. 

tempi  _d, 

templ_I_tr,  templ_I_tg,  templ_I_tb, 
tempi _I_sr,  templ_I_sg,  tempi _I_sb, 
ray_top  ); 

end; 

•*nd: 


{.PA} 

f  •*«■***»****«***»#***»*»**««******************»«*«**«»********************  j 
i  ************* ************ *  *  ***  at  m  9  VlAlN  ******  *  999  9  **  9  *«»**  x*»«x*x«x  xxx  *  *x  *  j 

begin 


{  INPUT  FILE  } 
assign  (sysin,  ’picture5.pas’); 
reset  (sysin); 


{  OUTPUT  FILE  } 
assign  (outfile,  ’picO.dta’); 
rewrite  (outfile); 


{  ************************************************************************** 

{  SET  UP  COUNTERS  } 
objcntr  =1; 
subobj  cntr  :=  1; 
light_cntr  :=  1; 

Dart  -ntr  -  l: 
ooiy  ntr  -  1: 
vert  ice  -ntr  —  l: 

{  CREATE  POINTERS  TO  RECORDS  } 
new(objcurr); 
new(subobjcurr); 
new(cpartcurr); 
new  (light  _current); 
new(poly_curr); 
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{  ***  READ  IN  DATA  FILE  ”*  } 


{  READ  IN  PICTURE  RECORD  } 
readln  (sysin,  picture. view  position  x); 
readln  (sysin,  picture. view  position  y); 
readln  (sysin.  picture. view  position  z); 
readln  Ssvsin.  picture. oackgrouna  roior  ri; 
readln  (sysin.  picture. background  color  g); 
readln  (sysin,  picture. background  color  b); 
readln  (sysin,  picture. screen_max  x); 
readln  (sysin,  picture.screen_max_y); 
readln  'svsin.  picture  ambient  r); 
readln  (sysin,  picture.ambient  g); 
readln  (sysin,  picture.ambient  b); 
readln  (sysin,  picture.no  zero); 
readln  (sysin,  picture. global_refraction_index); 
readln  (sysin,  picture. num  lights); 
picture. lights  :=  light  ^current; 

while  picture. num  lights  >  0  do  begin 

{  READ  IN  LIGHT  DATA  } 

readln  (sysin, light_current '(light  cntr  . I  rj; 
readln  (sysin, light_current  '[light  cntr  . I  gj; 
readln  (sysin, light  current  '[light  cntr. I  b); 
readln  (sysin, lightcurrent  *  (light  _cntrj .  light_x) ; 
readln  (sysin, light  current ' [light _cntrj . light _y) ; 
readln  (sysin, light  current  '(light  cntrj.light  z); 
readln  (sysin, light  current  '[light  cntr]. dimensionl); 
readln  (sysin, light  current  '[light_cntrj.dimension2); 
lightcntr  :=-  light  cntr  +  1; 
picture  num  lights  .=  picture. num  lights  -  1; 
end; 

readln  (sysin,  picture. num  objs); 
object  loop _cnt  :=  picture. num  objs; 
picture. objects  :=  obj  curr; 


while  object  loop  cnt  >  0  do  begin 

{  READ  IN  OBJECT  DATA  } 

readln  (sysin,  obj_curr '[obj _cntr|. opcode); 

readln  'svsin.  obj  curr'  obj  cntr  .obj  bsphere  radius): 

■eauln  svsin.  > o |  :urr  '  on)  cntr:. obj  bspnere  s); 

-eat!  In  svsin.  >oj  curr  '  ooj  -ntr  ooj  osohere  v(; 

readln  (sysin,  ooj_curr  iooj_cntrj.ooj  ospnere  i); 

readln  (sysin,  obj_curr '(objcntrj.num  sub  objects); 
subobj  _loop_cnt  :=  obj_curr  '[obj  cntrj.num  sub  objects; 
picture  objects  *[obj_cntr]. sub  objects  :=  subobj  curr; 

while  subobj _loop  cnt  >  0  do  begin 

{  READ  IN  SUBOBJECT  DATA  } 

readln  (sysin,  subobj_curr ' [subobj _cntr]. subobj  type); 
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readln  (sysin,  subobj_curr '!subobj_cntrj.sub_bsphere  radius); 
readln  (sysin,  subobj_curr'[subobj_cntr|.sub_bsphere  x); 
readln  (sysin,  subobj _curr '[subobj _cntr|. sub  bsphere  y); 
readln  (sysin,  subobj  curr  *jSubobj_cntrj.sub_bsphere  z); 
readln  (sysin,  subobj  curr  *[subobj_cntrj.num_common_parts); 
cpart  loop  cnt  :=  subobj  curr'fsubobj  _cntr|.num  common  parts; 
picture. objects  '  obj  cntr  .sub  jibjects  *  subobj  _:ntr  . common  parts  :  = 
cpart  _curr; 

while  cpart_loop_cnt  >  0  do  begin 

{  READ  IN  COMMON  PART  DATA  } 

readln(sysin,  cpart_curr  '?cpart_cntr!.K  _arl; 
readln(sysin,  cpart_curr  *[cpart_cntrj.K_ag); 
readln(sysin,  cpart_curr  *[cpart_cntrJ.K_ab); 
readln(sysin,  cpart_curr  '(cpart_cntr].K_dr); 
readln(sysin,  cpartcurr  ‘(cpart_cntr|.K_dg); 
readln(sysin,  cpart  curr  *[cpart_cntr].K_db); 
readln(sysin,  cpart  curr '[cpart_cntr].K_sr); 
readln(sysin,  cpart  curr '  cpart  cntri.K  sg); 
readln(sysin,  cpart_curr  '|cpart_cntrj.K  sb); 
readln(sysin,  cpart  curr* [cpart  cntr]. K  tr); 
readln(sysin,  cpart  curr  *[cpart_cntr).K_tg); 
readln(sysin,  cpart  curr  *[cpart_cntrj.K  tb); 
readln(sysin,  cpart_curr  '|Cpart_cntr;.obj  ret'raction_indexj; 
readln(sysin,  cpartcurr '[cpart_cntr].obj_phong_exp); 

{  CHECK  TO  SEE  IF  SUBOBJECT  IS  A  SPHERE  OR  A  POLYGONAL  OBJECT  } 
if  (subobj_curr*[subobj_cntr].subobj_type  =  1)  then  begin 

readln  (sysin ,  cpart  curr  *  [cpart  cntr] .  n  urn  poly  gons) ; 
poly  loop  cnt  :=  cpart_cuiT*[cpart_cntr].num_polygons; 
picture. objects  '[obj_cntr].sub_objects "[subobj  cntr]. 
commonparts '[cpartcntr]. polygons  :=  poly_curr; 

while  poly  loop  cnt  >  0  do  begin 
{  READ  IN  POLYGON  DATA  } 

readln(sysin,  poly  curr  *[poly_cntrJ.num_vertices); 
vertice  loop  cnt  :=  poly _curr '[poly  cntrj.num  vertices; 

while  vertice  loop  cnt  >  0  do  begin 
{  READ  IN  VERTICE  DATA  ] 

readlnlsvsin,  poly  curr'^poly  cntri  .vertice  Uvertice  cntri); 
reatllnlsysin.  ooiv  '■urr  poly  :ntri. vertice  y  vertice  ••  n cn ) ; 
readlnlsvsin.  joiv  'urr'  poiy  -ntr*. vertice  z  vertice  cntr  ); 

vertice  cntr  :=  vertice  cntr  -r  l; 
vertice_loop_cnt  :=  verticejoop  cnt  -  1; 
end; 

vertice  cntr  :=  1; 

readln(sysin,  poly  curr '[poly  cntr], surface  normal  x); 
readln(sysin,  poly  curr  "(poly  _cntr[. surface  normal_y); 
readln(sysin,  poly_curr '[poly _cntr|.surface_normal_z); 
poly  cntr  :=  poly  cntr  +  1; 
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poly_loop_cnt  :=  poly  _loop_cnt  -  1 

end; 

poly  cntr  :=  1; 
new(poly_curr); 

end;  {  IF  j 

cpart  cntr  :=  cpart  cntr  1; 
cpart  loop  cnt  ;=  cpart  loop  cnt  -  1 

end; 

cpart_cntr  :=  1: 

new  (cpart  _curr); 

subobj  cntr  :=  subobj  cntr  +  1; 

subobj  loop  cnt  :=  subobj  _loop_cnt  -  1 

end; 

new  (subobj  _curr); 

subobj  cntr  =  l; 

obj  cntr  :=  obj  cntr  -t-  1; 

object  loop  cnt  :=  object  loop  cnt  -  1 


{PA} 


{  SET  RAY  STACK  POINTER  } 

ray  top  :=  nil; 

pixel  z  :=  initial  pixelz; 

{  RASTER  SCAN  LOOP  } 

for  pixely  :=  I  to  picture.screen_max_y  do  begin 
for  pixel  x  :=  1  to  picture. screen  max  x  do  begin 

{  DETERMINE  VIEW  RAY  DIRECTION  } 
x  :=  pixel  x  -  picture.view_position_x; 
y  ;  =  pixel_y  -  picture. view_position_y; 
z  :=  pixel  z  -  picture. view  position  z; 

dist  sartfsar(x)  -  sqrfyl  -  sqr(z)); 

'ON VERT  IT  TO  V  JMT  VECTOR 
inux  —  x  list: 
unity  :=  y  /  dist; 
unitz  ;=  z  /  dist; 

{  INITIALIZE  VIEW  RAY  } 

initial_ray  _type  :=  view; 
initial_ray_origin_x  :=  pixe)_x; 
initial_ray_ongin_y  :=  pixel_y; 
initial  ray  origin  z  :=  pixel  z; 
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imtial_ray  _vector  x  :  =  unitx; 
initial_ray  vector  y  :=  unity: 
initial  ray  vector  z  :=  unitz; 
imtial_ray_stype  :=  none; 
initial  intersection  flag  :=  false; 
initial  obi  idx  :=  0: 

;mciai  -iuoobj  cix  :  =  0; 

initial  opart  idx  :=  0; 

initial_polygon_idx  :=  0; 

initial_intersection  x  :=  0.0; 

initial_intersection_y  :=  0.0; 

initial  intersection  z  —  0.0: 

initial_d  :=  0.0; 
initial_I_tr  :=  0.0; 
initial_I_tg  :=  0.0; 

initial  I _ tb  :=  0.0; 

initial_I_sr  :=  0.0; 
initial  I  sg  :=  0.0; 

'.nitiai  I  so  :=  0.0: 


& 


ray  generation  number  :=  0; 


push!  initial  ray  type, 

initial  ray _ongin_x, initial _ray  origin _y, initial  ray  origin  z. 

initial  ray  vector  x,initial_ray_vector  y, initial  ray  vector  z, 

initialraystype, 

initial  intersection  flag, 

initial_obj_idx, 

initial  subobj  idx, 

initial  cpart  idx, 

initial  polygon  idx, 

initial_intersection  x, 

initial  intersect ion_y, 

initial  intersection _z, 

initial_d, 

initial_I_tr,  initial _ I  tg,  initial  I  tb, 

initial_I_sr,  initial_I  sg,  initial  I_sb, 
ray  top  ); 


{  BEGIN  RAY  TRACING  LOOP  } 


a 


oop  'current  "av  rype. 

current_ray _ongin_x,current_ray  origin  y,current_ray  origin  z, 

current_ray  vector  x, current  ray  vector  y, current  ray  vector  z, 

current_ray  stype, 

current_intersection  flag, 

current _obj  idx, 

current_subobj  idx, 

current _cpart  idx, 

current _polygon  idx, 
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current  intersection  _x. 
current  intersection_y . 
current  intersection  z. 
current_d, 

current  I_tr,  current_I_tg,  current_I  tb, 
current  I  sr.  current  !  «g,  current  I  sb. 

ray  top  j: 

if  (current  intersection_flag)  then  begin 

{ 

IF  THIS  FLAG  IS  SET  THEN  THIS  RAY  HAS  ALREADY  BEEN  THROUGH  THE  RAY 

TRACING  PROCESS  AND  HAS  HIT  AN  OBJECT. 

} 

cpart  path  :=  picture. 

objects '[current_obj_idx]. 

sub  objects '[current  subobj  idx] . 

common_parts; 

subobj  path  —  picture. 

objects  '[current_obj_idxj. 
sub_objects; 

if  subobj  path  “  current  subobj  idx  .subobj  type  =  0  then  begin 

{ 

IF  THE  SUBOBJECT  TYPE  IS  A  SPHERE  THEN  THE  SURFACE  NORMAL  AT  THE  POINT 
OF  INTERSECTION  MUST  BE  CALCULATED  SINCE  IT  CAN  NOT  BE  STORED.  THE  SURFACE 
NORMAL  IS  DETERMINED  FOR  THE  PLANE  TANGENT  TO  THE  SPHERE  AT  THE 
INTERSECTION  POINT. 

} 

surface  ^normal  :=  current  intersectionx  - 

subobj  path  “ [current  subobj  idx] . 
sub_bsphere_x; 

surface  normal  y  :=  current  intersection  y  - 

subobj  path '[current  subobj_idx]. 
sub_bsphere_y; 

surface_normal_z  :=  current _intersection_z  - 

subobj_path  '[current  subobj_idx]. 
sub_bsphere_z; 

•list  -  '  sqrtl  sqrfsurface  normal  x)  - 
sqr(  surface  normal  yj  - 
sqrisurface  normal  zi)): 

{  THIS  RAY  IS  THEN  CONVERTED  INTO  A  UNIT  VECTOR  } 
surface  normal  x  :=  surface_normal  x  /  dist; 
surface  normal_y  :=  surface_normal_y  /  dist; 
surface_normal_z  :=  surface_normal  z  /  dist; 
end 


IF  THE  SUBOBJECT  IS  A  POLYGONAL  OBJECT  THEN  THE  SURFACE  NORMALS  FOR  EACH 
OF  THE  POLYGONS  OF  WHICH  IT  IS  COMPOSED  IS  RETRIEVED  FROM  IT’S  RECORD. 

} 

surface  _normal_x  :=  cpart_path  *[current_cpart_idx,. 

polygons ‘[current  polygon  idx). 
surface  normal  _x; 

surface_normal  y  :=  cpart_path  *jcurrent_cpart_idx  . 

polygons  ‘(current_polygon_idx|. 
surface_normal_y; 

surface_normal_z  :=  cpart_path ‘Icurrent  cpart_id.x'. 

polygons ‘[current_polygon  idx j . 
surface_normal_z; 

end; 


{ 

PROCEDURE  FOR  DETERMINING  THE  INTENSITY  OF  LIGHT  AT  EACH  INTERSECTION 
POINT. 

} 

12  (currentd, 

currentl  tr,  current_I_tg,  currentltb, 
current  I  sr,  current_I_sg,  current_I_tb, 
picture. num  lights, 

picture. ambient  r,  picture. ambient_g,  picture. ambient  b, 

cpart  path  ‘ [current  cpart  idx]. K_ar, 

cpart  path  ‘[current  cpart  idxj.K  ag, 

cpart  path  ‘ [current  cpart  idx]. K  ab, 

cpart  path  ‘[current  cpart  idxj.K  dr, 

cpart  path  *  [current  cpart  idxj.Kdg, 

cpart  path  ‘[current  cpart  idxj.K  db, 

cpart_path  ‘[current  cpart  idxj.K  sr, 

cpart_path  ‘[current_cpart  idxj.K_sg, 

cpart _path  ‘[current_cpart_idx].K_sb, 

cpart  path  ‘[current  cpart  idxj.K  tr, 

cpart  path ‘[current  cpart  idx].K_tg, 

cpart  path  ‘[current_cpart_idxj.K_tb, 

current  _ray_type, 

raytop, 

picture. lights, 

picture. objects, 

cpart  path  ‘'current  cpart  idx  obj  ohong  cxd. 

urrent  -av  .ector  \. 
current  rav  vector  y. 

current_ray_vector_z, 
current  intersection  x, 
current  intersection  y, 
current  _intersection_z, 

surface_normal  x,  surface  normal  y,  surface  normal  z, 
ray  generation_number, 
intensity  _red, 
intensity  green, 


intensity  blue); 
end 

else  begin 

{  CHECK  FOR  POSSIBLE  INTERSECTIONS  OF  CURRENT  RAY  WITH  OBJECTS  IN  SCENE  } 

check  for  tncerseccionjcurrenc  rav  vector  x, 
current  ray  vector  y, 
current  _ray_vector_z, 
current  _ray_origin_x, 
current_ray_origin_y, 
current  ray  origin  _z, 
picture. objects, 
current  intersection  x, 
current_intersection  y, 
current_intersection  z, 
current  _obj_idx, 
current  _subobj  _idx, 
current  cpart  idx, 
current_polygon_idx, 
current_intersection  flag); 


{  SET  UP  PATHNAMES  TO  USE  AS  SHORTHAND  } 

subobjpath  :=  picture. 

objects  '(current_obj_idxj. 
subobjects; 

cpart  path  :=  picture. 

objects  *  [current  _obj  idx  j . 
subobjects  '(current_subobj_idx|. 
commonparts; 

if  subobj  path '[current_subobj_idx].subobj_type  =  0  then  begin 
surface_normal_x  :=  current  intersection  x  - 

subobj_path '[current  subobj  idx). 
subbspherex; 

surface_normal_y  :=  current  intersection  y  * 

subobj_path '(current  subobj  idx]. 
sub_bsphere_y; 

surface  normal  z  :=  current  intersection  z  - 

subobj  path  "current  subobj  idx1. 
iiio  bsDherp  z: 

list  :  -  sort  I  sort  surface  normal  x)  - 
sqr;surface  normal  y)  -r 
sqr(surface_normal_z))); 

surface_normal_x  :=  surfacenormalx  /  dist; 
surface  normal_y  :=  surface  normal  y  /  dist; 
surface  normal_z  :=  surface  normal  z  /  dist; 
end 

else  begin 
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surface  normal  x  :=  cpart  path ''current  cpart  idx!. 

polygons '[current  polygon  idx]. 
surface  normal  x; 

surface_normal_y  :=  cpart  path  *[current_cpart_idx]. 

polygons  "^current  polygon  idxl. 
surface  normal  y; 

surface_normal_z  :=  cpart  path 'jcurrent_cpart_idx]. 

polygons  *  [current  _polygon_idx]. 
surface  _normal  _z ; 

end: 

if  (current_intersection_flag)  then  begin 

{  IF  THERE  HAS  BEEN  AN  INTERSECTION  THEN  CONTINUE  TRACING  THE  RAY  } 

if  (stack_exceeded(  ray_generation_number, 

maximum  size  of_stack))  then  begin 

(  IF  THE  STACK  IS  ALREADY  FULL  THEN  CALCULATE  INTENSITY  AT  LAST  NODE  \ 

16  (current  d, 

current  I_tr,  current  I  tg,  current_I  tb, 

current  l  sr,  current  I  sg,  current  _I_tb. 

piccure.num  lights, 

picture.  ambient_r, 

picture,  ambient  _g, 

picture,  ambient  _b, 

cpart  path  "(current  cpart _idx|.K_ar, 
cpart  path  "[current  cpart  idx]. K_ag, 
cpart  path  "(current  cpart  idx) .K  ab, 
cpartpath  *  [current  cpartidxj.Kdr, 
cpart  path  *  [current  _cpart_idx].K_dg, 
cpart_path  "[current  cpart  idx].K_db, 
cpart_path  "[current  cpart  idx].K_sr, 
cpart _path  '  [current  _cpart_idxj.K_sg, 
cpart _path  " [current _cpart_idx].K_sb, 
cpart_path  "(current_cpart  idxj.K_tr, 
cpart_path  "(current  cpart  idx].K_tg, 
cpart  path  “(current_cpart_idx].K_tb, 
current  ray_type, 
ray  _top, 
oicture. lights, 
oict.  are.  objects. 

part,  patn  current  -part  ;dx  oDj  phong  ^*xp. 

current_ray  _vector_x, 
current  ray_vector  y, 
current _ray_vector  z, 
c  urre  n  t  _in  tersec  t  ion  _x , 
current  in tersection  y, 
current  intersection  z, 
surface  normal  x, 
surface  normal  y, 


surface  normal_z. 
intensity  _red, 
intensity  _green, 
intensity  _blue) 
end 

else  begin 

'  IF  THERE  WAS  AN  INTERSECTION  AND  THE  STACK  WAS  NOT  FULL  THEN  CALCULATE 
THE  DISTANCE  BETWEEN  THE  RAY’S  ORIGIN  AND  POINT  OF  INTERSECTION  AND  PLACE 
THE  RAY  BACK  ON  THE  STACK. 

} 

current_d  :=  (sqrt(sqr(current_intersection_x  - 
current_ray_origin_x))  -r 
(sqr(current_intersection  x  - 
current_ray_origin  x) )  + 

(sqr(current  intersection  _x  - 
current  _ray  _origin  _x) ) ) ; 

push(  current_ray_type, 
current  ray  origin_x, 
current  ray_origin_y, 
current  ray  origin  _z. 
current  ray  vector  x, 
current  _ray_vector_y, 
current  _ray_vector_z, 
current  ray  sty pe, 
currentintersectionflag, 
current  obj  idx, 
current  subobj  idx, 
current  cpart  _idx, 
current_polygon  idx, 
current  _intersect  ion  _x, 
current  _intersection_y, 
current  intersection  z, 
currentd, 
current_I_tr, 
current_I_tg, 
current  I  tb, 
current  I  sr, 
current  _I_tg, 
current  I  fb. 

-ay  t.op  j; 

{  DETERMINE  IF  A  REFLECTED  RAY  WAS  CREATED  AND  IF  SO  CALCULATE  IT  } 
calculate  reflected  ray  (current  ray  vector  x, 
current_ray_vector  y, 
current_ray_vector  z, 
surface  normal  x, 
surface  normal  y, 
surface  normal  z, 
reflected  ray  x, 
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reflected  ray_y, 
reflected  ray_z, 
reflected  ray  ); 


f  DETERMINE  IF  A  REFRACTED  RAY  W  AS  CREATED  AND  IF  SO  CALCULATE  IT  ) 

calculate  ret'raccea_ray  (current  ray  vector  x. 

current  _ray_vector_y, 
current  _ray_vector_z, 
surface_normal_x, 
surface  normal  y, 
surface  normal  z, 
cpart_path “[current  cpart  idx]. 
obj_refraction_index, 
picture. global_refraction  index, 
refracted  ray  x, 
refracted_ray_y, 
refracted  ray  z, 
refracted  ray  i : 

raygenerationnumber  :=  raygenerationnumber  +  1; 
if  reflected  ray  then  begin 


IF  A  REFLECTED  RAY  WAS  CREATED  THEN  INITIALIZE  IT  AND  PUSH  IT  ON  THE 
STACK. 

} 

sourceraytype  :=  current  ray  type; 

dist  :=  (  sqrt(sqr(reflected_ray  x)  + 
sqr(reflectedrayy)  + 
sqr(reflected_ray_z))); 

{  CONVERT  REFLECTED  RAY  TO  A  UNIT  VECTOR  } 
unitx  :=  reflectedrayx  /  dist; 
unity  :=  reflected  ray_y  /  dist; 
unitz  :=  reflected _ray_z  /  dist; 

initialj  ay  type  —reflected; 
initialrayoriginx  :=  current_intersection  x; 

Initial  ray  origin  v  current  intersection  y: 
initial  rav  origin  z  :=  current  n terser t ion  :: 
initial  rav  vector  •<  :  -  unitx: 
initial_ray_yector  y  :=  unity; 
initial  ray  vector  z  :=  unitz; 
initial_ray_stype  :=  source  ray  type; 
initial_intersection_flag  :=  false; 
initial  obj  idx  :=  0; 

initial  subobj_idx  :=  0; 

initial_cpart_idx  :=  0; 

initial  polygon  idx  .=  0; 


initial  intersection _x  :=  0.0: 

initial  intersection _y  :=  0.0; 

initial  intersection _z  :=  0.0: 

initial_d  :=  0.0; 

initial_I_tr  :=  0.0; 

initial  I  _tg  :=  0.0: 

miciai  I  tb  :-=  0.0; 

initial  I  sr  :=  0.0; 

initial_I_sg  :=  0.0; 

initial_I_sb  :=  0.0; 

push(  initial_ray_type, 
initial_ray_origin_x, 
initial_ray_origin_y, 
initial  ray_origin  z, 
initial  ray_vector_x, 
initial_ray_vector_y, 
initial  ray  _vector_z, 
initial  _rav  stvpe. 
initial  intersection  _tiag, 
initial_obj_idx, 
initial  subobj_idx, 
initialcpart  idx, 
initial_polygon_idx, 
initialintersectionx, 
initialintersectiony, 
initialintersectionz, 
initiald, 

initialltr,  initialing,  initial  I  tb, 
initialise  initiallsg,  initiallsb, 
ray  top  ); 

end; 

if  refracted_ray  then  begin 

{ 

IF  A  REFRACTED  RAY  WAS  CREATED  THEN  INITIALIZE  IT  AND  PUSH  IT  ON  THE 
STACK. 

} 

source  rav  type  :  -  current  ray  fvpe; 

disc  !  iqrtlsorl  refracted  ray  xi  - 

sqr(refracted_ray_y)  + 
sqr(refracted  ray  z) ) ) ; 

{  CONVERT  IT  TO  A  UNIT  VECTOR  } 

unitx  :=  refracted _ray  x  /  dist; 
unity  :=  refracted_ray  y  /  dist; 
unit*  :=  refracted  ray  z  /  dist; 


inmai_ray  type  :=  refracted; 

initial _ray_origin_x  :=  current  intersection  x; 

initial_ray _origin_y  :=  current  intersection  y; 

inmal_ray  origin  z  :=  current  intersection  z; 

initial_ray_vector_x  :=  unitx; 

initial  ray  vector  y  :=  unity; 

initial _ray  vector  z  :  -  unitz; 

initial  _ray_stype  :=  source  _ray  type; 

initial_intersection_flag  :=  false; 

initialobjidx  :=  0; 

initial_subobj_idx  :=  0; 

initial  cpart  idx  .=  0; 

initialpolygonidx  :=  0; 

Lnitial_intersection  x  :=  0.0; 

initial_intersection_y  :=  0.0; 

initial_intersection_z  :=  0.0; 

initial_d  :=  0.0; 

initial_I_tr  :=  0.0; 

initialling  0.0; 

initiai_i_tt>  :=  0.0; 

initial_I_sr  :=  0.0; 

initial_I_sg  :=*  0.0; 

initial  I  so  :=  0.0; 


push(  initial  ray  type, 
initial_ray_origin_x, 
initialrayoriginy, 
initial  ray  origin _z, 
initial  ray _vector_x, 
initial  ray  vector  y, 
initialrayvectorz, 
initial  _ray_stype, 
initialintersection  flag, 
initialobjidx, 
initial  subobj  idx, 
initial  cpart  idx, 
initial_polygon_idx, 
initial  intersection  x, 
initial_intersection_y, 
initial  intersection  z, 
initial  _d , 

initial  1  tr.  initial  l  tg,  initial  1  tb. 
nitiai  1  'T.  nitiai  i  ->g,  initial  i  -o. 
-av  'on 


end; 

er  d; 
end 

else  begin 
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if  fcurrem  rs\ 


vie*  then  be<m 


'•V  [>* 

{ 

IF  THERE  WAS  NO  INTERSECTION  \ND  I  HE  (  I  RRENT  R  AV  IS  THE  VIEW  RAY  THEN  SET 
THE  OUTPUT  INTENSITY  TO  THE  BACKORul  NO  INTENSITY 

} 

intensity  r»d  •■=  picture  backereuind  <*-->l*>r  r 
.ntensuv  ,>reen  -  ucture  :><i<  K«r'M(n.l  ■  nor  -> 
intensity  blue  —  picture  baoitgroiind  coior  b 

end; 

end; 

end; 

until  (stack_empty(ray_top)); 

{  OUTPUT  THE  FINAL  INTENSITY  } 
display  pixel(  intensity  red, 
intensity  green, 
intensity  blue, 
ptxelx, 
pixel  _y, 
pixel  _z  ); 

end 

end: 

dose  (sysin); 
close  (output); 


end.  {  MAIN  } 


APPENDIX  B  -  INPUT  FILE 


100  ,  x/  view  position  {  PICTURE  } 

100  /y/ 

1000  /z/ 

0.0  /red/  background  light 
0.0  /green/ 

1.0  /blue,/ 

200  /x/  screen 

200  /y/ 

1.0  / red/  ambient  intensity 

1.0  /  green/ 

1.0  /blue/ 

1.0  /global  refraction  index/ 

1  /num  lights/ 

1.0  /red/  intensity  of  light  source 
1.0  /green/ 

1.0  /blue/ 

0.0  /x/  position  of  light  source 

20.0  /y/ 


0.0  /■/ 

0.0  /dimensionl/ 

0.0  /dimension2/ 

3  /numobjects/ 

9999  /opcode/  {  OBJECT  1  } 

35.0  /radius  of  object’s  bounding  sphere/ 

110.0  /x/ center  of  bounding  sphere 
0.0  /y/ 

-50.0  III 

1  /number  of  subobjects/ 

1  /subobject  type/  {SUBOBJECT  1} 

35.0  /radius  of  subobjects  bounding  sphere/ 

110.0  /x/  center  of  bounding  sphere 
0.0  /y/ 

-50.0  /z/ 

1  /num-common-parts/  {  COMMON  PART  1  } 


0.8 

/Ka-red/ 

ambient  coefficient 

0.0 

/Ka-green/ 

1.0 

Ka-olue 

').$ 

Kd-recl 

iitfuse  ■oerficient 

0.0 

,  Ka-green, 

0.0 

/Kd-blue/ 

0.8 

/Ks-red/ 

specular  coefficient 

0.8 

/Ks-green/ 

0.8 

/Ks-b!ue/ 

0.0 

/Kt-red  / 

transmission  coefficient 

0.0 

/Kt-green/ 

0.0 

/Kt-blue/ 
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/obj-refraction- index/ 
/obj-phong-speculir-exponent/ 
/num-polygons/ 

/num-vertices/  {  POLYGON  1  } 
/poly  1  ptl  / 


0.0 
200 
6 
4 

90.0 
20.0 
-30.0 

90.0  /polyl  pt2 / 

-20.0 
-30.0 

130.0  /polyl  pt3/ 

-20.0 
-30.0 

130.0  /polyl  pt4 / 

20.0 

-30.0 

0.0  /polyl  surface  normal/ 

0.0 
1.0 
4 

90.0 
-20.0 
-30.0 

90.0  /  poly 2  pt2 / 

-20.0 

-70.0 

130.0  / poly2  pt3/ 

-20.0 

-70.0 

130.0  /poly2  pt4/ 

-20.0 

-70.0 

0.0  /poly2  surface  normal/ 

-1.0 

0.0 

4  /num-vertices/  {POLYGON  3} 
130.0  /poly3  ptl/ 

20.0 


/num-vertices/  {  POLYOGN  2  } 
/poly2  ptl/ 


-30.0 

130.0  ,  poiy4  pt2/ 

20.0 
-TO  0 

>0.0  poiv4  pc3  - 
20.0 

-70.0 

90.0  /poly4  pt4/ 

20.0 

-30.0 

0.0  /poly4  surface  normal/ 

1.0 

0.0 

4  /num-vertices/  {  POLYGON  5  } 

130.0  /poly 5  pt  1  / 

20.0 

-30.0 

130.0  /poiyo  pt2 
-20.0 
-30.0 

130.0  poiy.i  pt3  / 

-20.0 

-70.0 

130.0  /poly 5  pt4/ 

20.0 

-70.0 

1.0  /poly5  surface  normal/ 

0.0 

0.0 

4  /num-vertices/  {  POLYGON  6  } 

90.0  /poly6  ptl/ 

20.0 

-70.0 

90.0  / poly6  pt2/ 

-20.0 

-70.0 

90.0  /poly6  ptS/ 

-20.0 

-30.0 

90.0  polv6  pt4 

3  0.0 

-30.0 

-1.0  /poiyo  surface  normal, 

0.0 

0.0 

9999  /opcode/  ***  OBJECT  2  *** 

175  /radius  of  the  objects  bounding  sphere/ 
100.0  /x/  center  point  of  the  bounding  sphere 

-100.0  /y/ 


1  num-subobjects 

1  'subobject-type '  {  SUBOBJECT  1  } 

175  radius  of  subobjects  bounding  sphere 
100.0  /x/  center  of  bounding  sphere 

-100.0  V 

-100  0  'z' 

L  numoer  common  parts,  i  COMMON  ? \RT  l  ; 

0.0  red  Ka  ambient  coefficient 
0.7  /green/ 

0.0  /blue/ 

0.0  / red/  Kd  diffuse  coefficient 

0.7  green 
0.0  /blue 

0.3  /red/  Ks  specular  coefficient 
0.8  /green./ 

0.8  /blue/ 

0.0  / red/  Kt  transmission  coefficient 

0.0  /green/ 

0.0  blue. 

0.0  /objects  refraction  index/ 

200  /Phong’s  specular  exponent/ 

1  /number  of  polygons/ 

4  number  of  vertices  {  POLYGON  1  } 

0.0  /  poiy  1  ptl, 

20.0 

•200.0 

0.0  /polyl  pt2/ 

0.0 

0.0 

200.0  /polyl  pt3 / 

0.0 

0.0 

200.0  /  poly  1  pt4/ 

20.0 

-200.0 

0.0  /polyl  surface  normal/ 

0.99 

0.1 

9999  /opcode/  ****  OBJECT  3  **** 

40  /radius  of  objects  bounding  sphere/ 

140  /x/  center  of  bounding  sphere 

30  y 

-  i  5u 

!  numoer  >t  suoooieets. 

0  /subobject  type/ 

40  /radius  of  subobjects  bounding  sphere/ 

140  /x/  center  of  bounding  sphere 

30  /  y/ 

-150  /*/ 

1  /number  of  common  parts/  {  COMMON  PART  1  } 
0.5  / red/  Ka  ambient  coefficient 

0.0  /green/ 


va 


/blue/ 

/red/  Kd  diffuse  coefficient 
/green/ 

/blue/ 

/red/  Ks  specular  coefficient 
/green  ' 
blue 

/red/  Kt  transmission  coeffii 
/green/ 

/blue/ 

/refraction  index  for  object/ 
/Phong's  specular  exponent 1 


specular  coefficient 


transmission  coefficient 


mrmmQ 
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